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Abstract 

We present exact calculations of the partition function of the g-state Potts model for general q and temper- 
ature on strips of the square lattice of width Ly = 3 vertices and arbitrary length with periodic longitudinal 
boundary conditions, of the following types: (i) {FBCy, PBCx) — cyclic, (ii) {FBCy,TPBCx) = Mobius, 
(iii) {PBCy,PBCx) = toroidal, and (iv) {PBCy,TPBCx) = Klein bottle, where FBC and {T)PBC refer to 
free and (twisted) periodic boundary conditions. Results for the Ly — 2 torus and Klein bottle strips are also 
included. In the infinite-length limit the thermodynamic properties are discussed and some general results 
are given for low-temperature behavior on strips of arbitrarily great width. We determine the submanifold 
in the space of q and temperature where the free energy is singular for these strips. Our calculations are 
also used to compute certain quantities of graph-theoretic interest. 
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1 Introduction 



The q-statc Potts model has served as a valuable model for the study of phase transitions and critical 
phenomena ||l|, |^. On a lattice, or, more generally, on a (connected) graph G, at temperature T, this model 
is defined by the partition function 

Z(G,g,z;)= ^e-/'^ (1.1) 

with the (zero-field) Hamiltonian 

n = -jY,^<y.<y, (1-2) 

where (Ji = 1, ...,q are the spin variables on each vertex (site) i £ G] (3 = {kBT)~^; and (ij) denotes pairs 
of adjacent vertices. The graph G = G{V,E) is defined by its vertex set V and its edge set E; we denote 
the number of vertices of G as n = n{G) — \V\ and the number of edges of G as e(G) = \E\. We use the 
notation 

K = f3J , a = w-i=e^, v = a-l (1.3) 

so that the physical ranges are (i) a > 1, i.e., v > corresponding to oo > T > for the Potts ferromagnet, 
and (ii) < a < 1, i.e., —l<v<0, corresponding to < T < oo for the Potts antiferromagnet. One defines 
the (reduced) free energy per site / = — where F is the actual free energy, via 

/({G}, g, v) = hm ln[Z(G, g, v)'/"] (1.4) 

n — >oc 

where we use the symbol {G} to denote lim„^oo G for a given family of graphs. In the present context, this 
n ^ oo limit corresponds to the limit of infinite length for a strip graph of the square lattice of fixed width 
and some prescribed boundary conditions. 

Let G' = (V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex set V and an edge 
set E' C E. Then Z{G, q, v) can be written as the sum [||-[|| 

Z(G,g,«) = q'^'^'^v^'^'^ 

G'CG 
n(G) e(G) 

= E Y.^rsq'-V^ (1.5) 

r=k{G) s=0 

where k{G') denotes the number of connected components of G' and Zrs > 0. Since we only consider 
connected graphs G, we have fc(G) = 1. The formula (^^) enables one to generalize q from Z+ to M_|_ 
(keeping v in its physical range). The formula ( |l.5| ) shows that Z{G,q,v) is a polynomial in q and v 
(equivalently, a) with maximum and minimum degrees indicated in eq. (|l.5| ). The Potts model partition 
function on a graph G is essentially equivalent to the Tutte polynomial and Whitney rank polynomial 

||], ijl^-Q for this graph, as discussed in the appendix. 

In this paper we shall present exact calculations of the q-state Potts model partition function Z{G, q, v) 
for arbitrary q and temperature for strips of the square lattice of width Lj, = 3 vertices and arbitrary length 
Lx, with periodic longitudinal boundary conditions of the following types: (i) {FBCy, PBCx) = cyclic, 
(ii) {FBCy,TPBCx) = Mobius, (iii) [PBCy^PBC^) = toroidal, and (iv) {PBCy,TPBCx) = Klein bottle, 
where FBC and [T)PBC refer to free and (twisted) periodic boundary conditions, respectively, and the 
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longitudinal and transverse directions are taken to be x (horizontal) and y (vertical) . Results for the Ly = 2 
torus and Klein bottle strips are also included for comparison. This is an extension of our previous exact 
calculations of Potts model partition functions for Ly — 2 strips of the square lattice ||l3|, |l^ and other 
lattices |15 including cyclic and Mobius boundary conditions. (Results for free longitudinal boundary 



conditions are in [|l4|-|7| and |j, yj.) 

There are several motivations for these exact calculations of Potts model partition functions for strips of 
the square lattice. Clearly, new exact calculations of Potts model partition functions are of value in their 
own right. From these, in the limit of infinite-length strips, one derives exact thermodynamic functions and 
can make rigorous comparisons of various properties among strips of different lattices. In particular, with 
these exact results one can study both the T — Q critical point of the g-state Potts ferromagnet and T = 
properties of the Potts antiferromagnet. Further, via the correspondence with the Tuttc polynomial, our 
calculations yield several quantities of relevance to mathematical graph theory. 

Various special cases of the Potts model partition function are of interest. One special case is the zero- 
temperature limit of the Potts antiferromagnet. For sufficiently large g, on a given lattice or graph G, this 
exhibits nonzero ground state entropy (without frustration). This is of interest as an exception to the third 
law of thermodynamics |2^, This is equivalent to a ground state degeneracy per site (vertex), W > 1, 
since Sq = fcslnH^. The T — Q (i.e., v = —1) partition function of the above-mentioned g-state Potts 
antiferromagnet on G satisfies 

Z{G,q,-l)^P{G,q) (1.6) 

where P{G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring the vertices 
of the graph G with q colors such that no two adjacent vertices have the same color [|| |l^, The 
minimum number of colors necessary for this coloring is the chromatic number of G, denoted xi^)- We have 

W{{G},q) = lim P{G,qf''' . (1.7) 

n — *oo 

A subtlety in the definition due to noncommutativity limits 

lim lim P(G,g)^/" 7^ lim limP(G,g)^/" (1.8) 

q — *qs n — >oo tl — ^oo q — ^qs 

at certain points qs was discussed in . 

As derived in , a general form for the Potts model partition function for the strip graphs Gs considered 
here, or more generally, for recursively defined families of graphs comprised of m repeated subunits (e.g. the 
columns of squares of height Ly vertices that are repeated L^ = m times to form an L^ x Ly strip of a 
regular lattice with some specified boundary conditions), is 

Z{{G,Uq,v)= cz,G.A^z,G.A<l^v)r (1-9) 

where Nz^g,,x, cz,Gs.jy ^-^^ ^z.Gs.j depend on the type of recursive family Gs (lattice structure and boundary 
conditions) but not on its length m. For special case of the T — antiferromagnet, the partition function, 
or equivalently, the chromatic polynomial P{{Gs)rm q) has the corresponding form 

P((G,)™,g)= cp,G..,(Ap,G.,,(g))™ . (1.10) 
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For sufficiently large integer q values, the coefficients in tfie cfiromatic polynomial ( 1.10| ) can be interpreted 



as multiplicities of eigenvalues of coloring matrices Tp |49[ , and the sums of these coefficients are thus sums 
of dimensions of invariant subspaces of these matrices. An analogous statement applies to the coefficients in 



the Potts model partition function (1.9), in terms of a generalized coloring matrix Tz- For a given family 



Gs of strip graphs wc shall denote these sums as 

Cz,G. = J2 "Z.G.,, (1-11) 

and 

Cp,G. = J2 ''P'G,,j ■ (1.12) 
J = l 

Below we shall give results for the determinants of the coloring matrices; following the notation in [55|] we 
have 

detTz(G.)= n (Az,G,,jr^'^=- (1-13) 

and 

detrp(G,)= H {Xp^G.^jT^-"-' ■ (1.14) 
i=i 

Let us define the shorthand notation 

Dp(G,) = detrp(G,) , DziGs) = detrz(G,) . (1.15) 

Here we distinguish the quantities Nz,Gs,\, -^p.Gs.a; cz,Gs,jj cp^G^.j! Az,GsJ: ^^^d Xp^Cs.j- Below, for 
brevity of notation, wc shall sometimes suppress the Z or P subscripts and the s subscript on Gs when the 
meaning is clear from context and no confusion will result. 

Using the formula ( |1.5| ) for Z(Gs, 9, v), one can generalize q from Z+ not just to but to C and a from 
its physical ferromagnetic and antiferromagnetic ranges 1 < a < oo and 0<a<ltoaeC. A subset of 
the zeros of Z in the two-complex dimensional space defined by the pair of variables (g, a) can form an 
accumulation set in the n — s- cx) limit, denoted B, which is the continuous locus of points where the free 
energy is nonanalytic. This locus occurs where there is a non-analytic switching between different \z,Gs,3 
of maximal magnitude and is thus determined as the solution to the equation of degeneracy in magnitude 



of these maximal or dominant Az,G,,i's 13, For a given value of a, one can consider this locus in the q 
plane, and we denote it as Bq{{G}, a). In the special case = (i.e., v = —1) where the partition function 
is equal to the chromatic polynomial, the zeros in q are the chromatic zeros, and i3g({G},a — 0) is their 
continuous accumulation set in the n — > c» limit. In previous papers we have given exact calculations of the 
chromatic polynomials and nonanalytic loci Bq for various families of graphs. Some early related papers are 
[^-[^. Our methods for calculating these chromatic polynomials were discussed in ||2^,j3^-^^ and some 
relevant bounds were given in p7[,||3^-p^. The families of of graphs of primary interest here are strips of 
regular lattices with periodic longitudinal boundary conditions; previous related papers on such families of 
graphs are |25|, [^-[^. With the exact Potts partition function for arbitrary temperature, one can 
study Bq for a 7^ and, for a given value of q, one can study the continuous accumulation set of the zeros 
of Z{G,q,v) in the a plane; this will be denoted Ba{{G},q). It will often be convenient to consider the 
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equivalent locus in the u = l/a plane, namely Bu{{G},q). We shall sometimes write Bq{{G},a) simply as 
Bq or B when {G} and a are clear from the context, and similarly with Ba and Bu- One gains a unified 
understanding of the separate loci Bq{{G}) for fixed a and Ba{{G}) for fixed q by relating these as different 
slices of the locus B in the space defined by {q, a) as we shall do here. As discussed earlier |13, 1^, in the 



infinite-length limit where the locus B is defined, for a given width and transverse boundary condition (free 
or periodic) B depends on the longitudinal boundary conditions but is independent of whether they involve 
orientation reversal or not. Thus, in the present context, for a given Ly, B is the same for the cyclic and 
Mobius strips, and separately for the torus and Klein bottle strips. 

Following the notation in and our other earlier works on Bq{{G}) for a = 0, we denote the maximal 
region in the complex q plane to which one can analytically continue the function T/F({G}, q) from physical 
values where there is nonzero ground state entropy as i?i . The maximal value of q where Bq intersects the 
(positive) real axis was labelled qc{{G}). Thus, region i?i includes the positive real axis for q > qc{{G}). 
Correspondingly, in our works on complex-temperature properties of spin models [|7^ we had labelled the 
complex-temperature extension (CTE) of the physical paramagnetic phase as (CTE)PM, which will simply 
be denoted PM here, the extension being understood, and similarly with ferromagnetic (FM) and antifer- 
romagnetic (AFM); other complex-temperature phases, having no overlap with any physical phase, were 
denoted Oj (for "other" ) , with j indexing the particular phase j?^, |75) . Here we shall continue to use this 
notation for the respective slices of B in the q and a oi u planes. Various special cases of Z{G, q, v) applicable 
for arbitrary graphs G were given in fl^ . Just as was true for the zero-temperature antiferromagnet, for the 
full partition function, at certain special points qs (typically qs = 0, l...,x(G)) one has 

lim lim Z(G,g,w)i/" 7^ lim lim Z(G, (j, u)^/" . (1.16) 

n — >oo q — >qs q — >qs n — >oo 



Because of this noncommutativity, the formal definition (1^) is, in general, insufficient to define the free 
energy / at these special points qs] it is necessary to specify the order of the limits that one uses in eq. (1.1£) 

The Potts ferromagnet has a zero-temperature phase transition in the ^ oo limit of the strip graphs 
considered here, and this has the consequence that for cyclic and Mobius boundary conditions, B passes 
through the T = point u = 0. It follows that B is noncompact in the a plane. Hence, it is usually more 
convenient to study the slice of B in the u = 1/a plane rather than the a plane. For the ferromagnet, since 
a — > oo as T — > and Z diverges like this limit, we shall use the reduced partition function Zr 

defined by 

Zr{{G,)rn, <Z, v) = O'^^^^ Z( (G, )„ , q, v) = U«''' Z {{G s) m , q, v) (1.17) 

which has the finite limit ^ q as T — s- 0. For a general strip graph {Gs)m of type Gg and length = m, 
we can write 

Zr{{GsU,q,v) = J2 cz,g.A^z,g.,jY^ ^ J2 cz.G..,(Az.G.,,.«)" (1.18) 

with 

Az,G„,,«=w^"'^=^'"^/"Az,G.,, . (1.19) 

For the Potts model partition functions of the width Ly strips of the square lattice with (i) cyclic or Mobius 
and (ii) torus or Klein bottle boundary conditions, the prefactor in eq. (1.19) is (i) u^^""^, and (ii) u^^", 
respectively. 
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For the following, it will be convenient to define some general functions. First, we define the following 
polynomial: 



P{Ck,q) 



k-2 



and P{Cn, q) is the chromatic polynomial for the circuit (cyclic) graph C„ with n vertices, given by 

P(C„, g) = iq- ir + (9 - 1)(-1)" . (1.21) 
Second, we define polynomials of degree d in g which are related to Chebyshev polynomials: 



M) - 



4) 



(1.22) 



where Un{x) is the Chebyshev polynomial of the second kind, defined by 

C/„(x) = ^(-l)^ r . ^ (2x)"-2^- (1.23) 



where in eq. ( 1.23| ) the notation \^ in the upper limit on the summand means the integral part of ^. These 
polynomials enter as coefficients in (1.9) and ( 1.10| ) for the cyclic and Mobius strips under study here. The 
first few of these coefficients are 

c(°) = 1 , c(i) = 9 - 1 , (1.24) 
c(2) ^ g2 _3^^ -^^ (^ 25) 

and 

c(3) = g3 5g2 _^ 6g - 1 . (1.26) 
The "falling factorial" used in combinatorics is defined as 

n-l 

(with 9(0) = 1), so that = g, q(2) = q{q - 1), 9(3) = q{q - 1)(9 - 2), etc. 

On each transverse slice of the width Ly cyclic and Mobius strips for Ly > 3, the (Ly — 2) interior 
vertices have degree (=coordination number) 4, while the two exterior vertices have degree 3. Let us define 
an effective degree or coordination number Ae/ j as 



'2\E 
n 

Then 



= hm (^) . (1.28) 



Ae//=4(l-^) ior sq,Ly,FBCy,iT)PBC^ . (1.29) 

Thus, for Ly = 3, ^eff = 10/3. The strips of the square lattice with torus and Klein bottle boundary 
conditions are A-regular graph with A = 4. (Here, a A-regular graph is defined as one in which each vertex 
has the same degree. A.) We proceed to describe our results. 
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2 Potts Model Partition Function for Cyclic and Mobius Ly = 3 
Strips of the Square Lattice 

2.1 General Structural Properties 

For cyclic and Mobius strips of the square lattice of fixed width Ly and arbitrary length Lx , the coefficients 
cz,Gs-j cp^Gs,j 6qs. (1.9) and (1.10) depend only on Ly and the boundary conditions. These coefficients 
are polynomials in g. In the case of cyclic boundary conditions, the coefficients of degree d in g have a unique 
form, namely S-'^^ given in ( 1.22| ), while in the case of Mobius boundary conditions, the coefficients of degree 
d are ic^''^. For the cyclic strip, let us denote the number of terms Xz,Gs,j ^-nd Xp,Gs,j in ^QS. (l.£) and 
( 1.10 ) respectively that have for their coefficient c^'^\ nz{Ly,d) and np{Ly,d). For the Mobius strip, we 
denote the number of terms Xz,Gs,j ^'^id Xp^G^j in eqs. ( pTol ) and (I.IC) that have for their coefficient ±c'-''\ 
nz{Ly, d, ±) and np{Ly, d, ±). In we proved a number of general theorems on the structure of Z{Gs, q, v) 
and P{Gs,q) that determined these numbers, as well as the overall number of terms Nz.Gs-X ^-^d Np^c^^.x- 
For our present purposes, we can apply the general formula 

^^^^.A ^ (2.1.1) 



Ly 



with Lj, = 3 to conclude that 



iVz.3,A = 20 . (2.1.2) 

Further, [psf 

(2d+l) ( 2Ly 



for < d < Ly, with nz{Ly, d) — for d > Ly] hence, for Ly = 3, it follows that 

nziS, 0) = 5 , nz(3, 1) = 9 , nz(3, 2) = 5 , nz{3, 3) = 1 . (2.1.4) 

For the Mobius strip, the analysis in gives, for the special case Ly — 3 

nz(3,0,+) = 6 ,nz(3,0,-) = 4 ,nz(3,l,+) - 6 ,nz(3,l,-) = 3 ,nz(3,2,+) = 1 (2.1.5) 

with other nz{3,d,±) ~ 0. 

In the special case of the zero-temperature Potts antiferromagnet, v = —1, from the general result p3] 



(Ly-j) 

one has, for Ly = 3 



Np-i,x = 10 . (2.1.7) 
Further, for the Ly — 3 cyclic and Mobius strips, 

np(3,0) = 2 ,7ip(3,l) = 4 ,np(3,2) = 3 ,np(3,3) = 1 (2.1.8) 

np(3, 0, +) = 3 , np(3, 0, -) = 2 , np(3, 1, +) = 3 , np(3, 1, -) = 1 , np(3, 2, +) = 1 (2.1.9) 
with other np(3, d, ±) — 0. 
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For the cyclic strip, ( |l.ll| ) and (|]T|) are ||5| 

Ly 



i=0 



and 



Cp,L, = E np{Ly,d)c^^^ = P{Tl^ , q) = q(q - l)^""! 



(2.1.10) 



(2.1.11) 



where r„ denotes the path or tree graph with n vertices. Hence, for the Ly — 3 cychc strip under study 
here, 

Cz,3 = g', Cp^s^qiq-l)^ . (2.1.12) 

The corresponding sums of coefhcients for the Mobius strips are, for the case of odd Ly relevant here, 

Cz^L„,A/6 = 9^^" + '^/" (2.1.13) 

and 



C 



P,Ly,Mb 



PiT{Ly+i)/2,q) = qiq- 1)" 



so that for Ly — 3, 



Cz,3,Mb ~ , Cp^s^Mb — q{q — i) 



(2.1.14) 
(2.1.15) 



2.2 Results for the Potts Model Partition Function 



Using a systematic iterative application of the deletion-contraction theorem and transfer matrix methods, 
we have calculated the exact Potts model partition function Z{G, q, v) for general q and v for the the cyclic 
and Mobius strips of the square lattice of width Ly = 3 and arbitrary length Lx — m. For quantities that are 
independent of L^, we shall use the abbreviations s3c and s3Mb for these strips. We find for these partition 
functions the results 



Z{sq,3 X m,FBCy,PBCx,q,v) = ^Cz^s3cji\z,s3cjiq,v)y 



and 



Z(sg,3 X m,FBCy,TPBCx,q,v) =J2^Z''3^^''ji\z,s3c,jiq:V))"' 
where (ordering the Xz,s3c.j^s by decreasing degree of the associated coefficients cz^sScj) 

Az,s3c,2 = v'^iq + v) 
■ + v{v + 5)± (^[q + v{v + 5)]2 - 4v{v + q){v + l)j 



A 



Z,s3c,(3,4) 



1/2- 



and 



A 



Z,s3c,{5,6) 



q + v{v + 3)± (^[q + v{v + 3)]^ - Av{v + q){v + 1)^ 



1/2- 



(2.2.1) 
(2.2.2) 

(2.2.3) 
(2.2.4) 
(2.2.5) 

(2.2.6) 
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Note that 



Az,s3c,3Az,s3c,4 = Az,s3c,5 Az,s3c,6 = "^(u + l)(w + q) 



The terms A2^s3c.j for j = 7, 8, 9 are roots of the cubic equation 



e + hie + /32e + /; 



33 



where 



and 



Further, 



/32 



hi^-v{2v + q){Zv + q + v'^) 
v^{v + q){v + l)(6w^ + Aqv + q'^ + v^) 

hs^-v'^v + qfiv + lf . 

Az,s346 ^ {v + q){v + . 



(2.2.7) 

(2.2.8) 

(2.2.9) 
(2.2.10) 

(2.2.11) 

(2.2.12) 



The terms Xz,s3c,j for 10 < j < 15 and 17 < j < 20 are roots of algebraic equations of degree 6 and 4, 
respectively, given as ( |9.1.l| ) and ( 9.1.8 ) in the appendix. The roots of the quartic equation ( 9.1.8 ) comprise 
the Xz,G,j for the 3 x strip of the square lattice with {FBCy, FBCx)- 

For the cyclic strip the coefficients in eq. ( [2.2.1 ) are, in terms of the c^''^ given above in eq. ( 1.22 ), 



and 



CZ,s3c..l = c^^^ 

cz,s3c.j = c'^^^ for 2 < j < 6 
cz.ssc.j = c^^^ for 7 < j < 15 

cz,s3c,j = c(") for 16 < j < 20 



„(2) 



For the Mobius strip, the coefficients in ( 2.2.2 ) are 

Cz,s3MbS — C 

cz,s3MKj - c(") for J = 5, 6, 17, 18, 19, 20 
cz,s3MKj =-c(°) for j = 2,3,4,16 
cz,s3MKj = c(i^ for 10 < j < 15 
cz.s3Mb.j = -c^^^ for 7 < j < 9 . 



(2.2.13) 

(2.2.14) 
(2.2.15) 

(2.2.16) 

(2.2.17) 

(2.2.18) 
(2.2.19) 
(2.2.20) 
(2.2.21) 



One easily checks that these results agree with the general structural formulas given above. 

From our calculations of Z{G, q, v) for both the Ly = 2 and Ly = 3 cyclic and Mobius strips of the square 
lattice, as well as the elementary calculation of the cyclic Ly = 1 case, we infer that the first two terms are 
consistent with the following general formula for arbitrary Ly and G = sq, Ly x L^, FBCy, {T)PBCx- 



A 



Z,G,1 



(2.2.22) 
(2.2.23) 
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with coefficients 

Cz,G.„..,i = c^'^''^ (2-2.24) 

and 

cz,G.„...2 = c(^''"'^ (2.2.25) 

for the cychc strips; for the Mobius strips, the corresponding coefficients can be determined by our general 
formulas in p5| . 

From these results we calculate the determinant 

20 

Dz{sq,3x L,,FBCy,PBC^) - l[{Xz,s3c,jr^-'''-' 

- (y- 1)39' 2/29^:^39' (2.2.26) 

where the variables x and y are defined in ( |9.2.3| ) and ( [j.2.4 ) in the appendix and are the variables of the 
Tutte polynomial. Comparing this with Dz{G) for the Ly ^ 1 and Ly ~ 2 strips we note that these 

cases Lj, = 1, 2, 3, are all described by the following formula 

Dzisq,Ly X L^,FBCy,PBC,) = (y - 1)^««"'' ^(^"-i)'"''"^^"'"''"' . (2.2.27) 

We also find, for the Mobius strip, 

20 

Dz(sq,3x L,,FBCy,TPBC^) = l[{Xz,s3c.jy^-''"-' 

= v39(9-l)(„^.i)2g(^_^^)3g 

= {y-lf''\^'^x^'^ . (2.2.28) 

Comparing this with Dz{sq, 2 x L^, FBCy, TPBC^) @, ||, we observe that both of these cases Ly — 
are described by the formula 

Dzisq,Ly X L^,FBCy,PBCo,) = (y - l)^''«^''"'y(^''-i)«x^''«^''"' . (2.2.29) 
2.3 V = —1 Special Case 

For the special case v — —1, i.e., the zero-temperature Potts antiferromagnet, where the partition fimction 
reduces to the chromatic polynomial, ten of the twenty Xz,s3c,j^s vanish and the others reduce to those 
for the Ly — 3 cyclic/ Mobius strip given in To see which \z,s3c,j^s vanish and which reduce to 

corresponding Xp.sSc.j^s, we recall the expressions for the chromatic polynomials: 

10 

P(s(?,3 X m,FBCy,PBC^,q)Y,cp,s3c.j{Xp,s3c,jr (2.3.1) 

10 

P{sq,3x m,FBCy,TPBC„q) - ^ cp,,3M6j (Ap,,3c,j)" • (2.3.2) 
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We order the terms Xp^sScj by decreasing degree of the associated coefficients cp^sScj for the cychc strip, as 
polynomials in q (this differs from the order and hence labelling used in |5^ ) . These terms are 



Ap,s3c,l — —1 

Ap,s3c,2 = 9 — 1 



P,s3c,3 



P,s3c,4 



and 



A 



P,s3c,5 



q~4 
q~2 



The terms Xp^sScj, j — 6, 7, 8, are the roots of the cubic equation 



with 



r + &3ir + ^32? + ^33 = 



631 = 2q^ -9(7+12 
632 = - lOq^ + 36q^ - 56q + 31 



(2.3.3) 

(2.3.4) 
(2.3.5) 
(2.3.6) 

(2.3.7) 

(2.3.8) 

(2.3.9) 
(2.3.10) 



Finally, we have 



A 



1 



P,s3c,(9,10) 



633 = -(g - l){q^ - 9q^ + 29q^ - 40<? + 22) . 



(q - 2)(g2 - 3g + 5) ± {((j2 - 5g + 7){q^ - 5q^ + Uq^ - 12q + 8)} 



1/2 



Thus, the correspondences are as follows. For v ^ —1, 

Az,s3c,l Xps3c,l 

and 

Az,s3c,2 Ap^s3c^2 ■ 

For the two pairs of quadratic roots, we have 

Az,s3c,3 A 



P,s3c,3 



A 



and 



Z,s3cA 



^Z,s3c,5 







\P,s3c,4 



A 



Z,s3c,6 







(2.3.11) 
(2.3.12) 

(2.3.13) 
(2.3.14) 



(2.3.15) 



(2.3.16) 



For V — —1 the cubic equation (2.2.8) reduces to ^^[^ + (9 — 2)-^] — 0, so, with the indicated labelling of 
roots, 



A 



Z,s3c,7 



A 



P,s3c,5 
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A 



Z,sSc,j 



for j ^8,9 . 



(2.3.17) 



The sixth order equation (|9.1.1 ) reduces to times the cubic equation ( 2.3.8 ) so that, with the indicated 
labelhng, 



A 



Z,s3c,j 



A 



Z.s3c,j 



Ap,s3c,j-4 for 10 < j < 12 



for 13 < j < 16 . 



(2.3.18) 



FinaUy, the quartic equation (9.1.S) becomes times a quadratic, yielding the reduction 

Az,s3cj Ap_s3cj-8 for j = 17, 18 

Az,.3c,j ^ for j = 19, 20 . 



(2.3.19) 



The sums and products of the eigenvalues of the coloring matrix for the Ly — 3 chromatic polynomials of 
the cyclic and Mobius strips were given in Note that, since some of the eigenvalues of the coloring 

matrix for the full Potts model vanish when v = — 1, one must extract these from the expressions ( ^.2.26| ) 



and ( 2.2.28 ) in order to recover the products for the corresponding chromatic polynomials. For example, for 
the cyclic Ly = 2> case, one has Q 



10 



J = l 



(2.3.20) 



2.4 V = {) Special Case 

For infinite temperature, u = 0, we have 

Az,s3c,i7 = for V = 

^z.sic.j = for V = Q and j ^ 17 
so that the partition function reduces to Z{sq,'i x m,FBCy, {T)PBCx,q,v = 0) = g". This can be seen 



(2.4.1) 
(2.4.2) 



directly from eqs. (1.1) and (1.2) or from eq. (1.5) 



3 Potts Model Partition Function for Ly = 2 and Ly = 3 Strips 
of the Square Lattice with Torus and Klein Bottle Boundary 
Conditions 

3.1 General Structural Properties 



For the strip with torus boundary conditions, (1.11) and (1.12) are |5£ 



(3.1.1) 
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as m ( ^.l.lOl ), and 

Cp^L=PiCr,q) (3.1.2) 



where C„ is the circuit graph and P{Cn,q) was given above in (1.21). Thus, for the Ly — 2 and Ly = 3 
strips with toroidal boundary conditions under study here, 

Cz.2^<l\ Cp,2 = qi2) = q{q - I) (3.1.3) 

and 

Cz,3 = g', Cp,3 = g(3)=9('?-l)(9-2) . (3.1.4) 

The corresponding sums of coefficients for the width Ly strips with Klein bottle boundary conditions (denoted 
k) are 

Cz,L,,k = gl^l (3.1.5) 



where [v] denotes the integral part of v, and 1 55 



Cp,L„k=0 (3.1.6) 

so that for Ly — 2, 3, 

Cz,2,k = q , Cz,3,k = g^ (3.1.7) 

and 

Cp,2,fc = Cp,3,A; = . (3.1.8) 

3.2 Results for the Potts Model Partition Function 

By the same methods as for the cyclic and Mobius strips, we have calculated the exact Potts model partition 
function Z(Gs , q, v) for general q and v for the the strips of the square lattice of width Ly = 3, arbitrary length 
Lx = TO, and torus or Klein bottle boundary conditions. For comparison, we also include our calculation of 
the partition function for the corresponding width Ly = 2 strips. For quantities that are independent of Lx, 
we shall use the abbreviations s2i, s2k, s3t, and s3k for these strips. 

3.3 Ly = 2 

Strictly speaking, the Ly ~ 2 strips of the square lattice with torus or Klein bottle boundary conditions 
are not proper graphs but multigraphs, i.e., they involve multiple (double) edges connecting vertical pairs of 
vertices. The number of vertices and edges are given by n = L^Ly and e — 4,Lx. The graphs in the family 
are A-regular with A = 4. We find for the partition function the results 

Nz,s2t,x = Nz,s2k,x = 6 (3.3.1) 

6 

Z{sq, 2 X TO, PBCy,PBC,,q, w) = ^ CZ,s2t,j{^Z,s2t,j{q,v))"^ (3.3.2) 

and 

6 

Z{sq,2xm,PBCy,TPBC,,q,v)=Y,cz.s2k.j{Xz.s2t.,3{q,v)r (3.3.3) 
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where 



A 



Z,s2t,(3,4) 



Az,«2t,i = (3.3.4) 
Az,s2f,2 = v{v + q) (3.3.5) 
^ g + 6u + + ± ^/R234 



(3.3.6) 



where 



where 



i?234 = 32^2 + 52.^ 8gw + 40v^ + 8v^ - 2qv^ + 24v^ (3.3.7) 
Az,s2t,(5,6) = \ [gv^ + ^qv + g2 + + qv^ +v^± ^/B^& ] (3.3.8) 



-R256 = (^^^ + + + ^qv^ + Qqv + q^){v'^ + + Av'^ - qv^ + 2qv + q^) . (3.3.9) 

We note that \z,sit,s^z,sit,i = v^{v + l)'^{v + q) while \z,sit,b^z,sit,Q = v'^{v + l)'^{q+v)'^. The corresponding 
coefficients for the torus and Klein bottle strips are 

cz,s2t,i = c^'^ (3.3.10) 

cz,«2fj = c(i) for .7-2,3,4 (3.3.11) 

cz,s2t,j = 1 for i = 5,6 (3.3.12) 

cz,s2k,i = -1 (3.3.13) 

Cz,s2k,2 = -c^'^ (3.3.14) 

cz,s2fc,,- =c« for j = 3,4 (3.3.15) 

cz,s2k,j =1 for j = 5,6 . (3.3.16) 
Prom these results we calculate the determinants 

6 

Dz{sq,2xL,,PBCy,PBC,) = n(^Wi)"'' 

i=i 

= u2«(«-l)(z; + l)29(i; + g)29 

= {y-lf^yxf' (3.3.17) 



and 



6 

.\Cz,s2k,j 



Dz{sq,2x L^,PBCy,TPBC^) = JJ(A 

= v^^''-'\v + lf\v + qf 

= [(y-l)2/]2V. (3.3.18) 

In the special case of the zero-temperature Potts antiferromagnet, v = —1, we have Np^g2t,\ = -^p,s2A;,a = 
4 and the corresponding reductions 

Az,.2t,i ^ 1 (3.3.19) 
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Az,s2t,2 1 — q 
Az,s2t,5 D4 = q"^ - 3q + 3 



and 



(3.3.20) 
(3.3.21) 
(3.3.22) 

(3.3.23) 



Az,s2t,j -> for j = 4,6 . 

In this limit the resuhant partition functions (chromatic polynomials) for the cases of torus and Klein bottle 
boundary conditions are identical to those for cyclic and Mobius boundary conditions, respectively. This is 
a special case of a general result that the chromatic polynomial for a graph G and a graph H which differs 
from G in having multiple edges between given sets of vertices are equal. 



3.4 L„ 



We find 



Recall that for the special case v 



Nz,s3t,x = 20 
Nz,s3k,X = 12 . 
-1 of the zero-temperature Potts antiferromagnet, 

5 . 



P,s3k,X — 



We calculate 



and 



20 



Z{sq, 3 X m, 

PBCy,PBC^,q,v)=Yl 



Z{sq, 3 X m, PBCy,TPBC^,q, v) = Y, cz.s3k.j (Az^sSfej (?, «))" 
where (ordering the Xz,s3t,j's by decreasing degree of the associated coefficients cz,s3t,j) 



A 



Z,s3t,j — Xz,s3c,j 



for 1 < j < 6 



A 



Z,s3t,{7,t 



y g + v{v + 6) ± (^{q + v{v + 6) f - 4:v{v + q){v + l)J 
Az,s3t,9 =v'^{v + l) 

Xz,s3t,w = w^(w + q){v + 1) . 
The next three terms coincide with three for the cyclic/Mobius strip: 

Az,s3tj = Az,s3cj-4 for j = 11, 12, 13 . 

The Xz,s3t,j for 14 < < 17 are roots of the quartic 

C'* + .941^^ + 342^^ + 343^ + 344 = 



1/2' 



(3.4.1) 
(3.4.2) 

(3.4.3) 
(3.4.4) 

(3.4.5) 
(3.4.6) 

(3.4.7) 

(3.4.8) 

(3.4.9) 
(3.4.10) 

(3.4.11) 
(3.4.12) 
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where 



= -v{25v'^ + 17v^ + v'^q + 8qv + + 6w* + q"^) 



(3.4.13) 



c,42 = v^{v + l)iq^ + iSv^ + v^q + A2v'^q + 2v^ + A2v'^ + 2,Qqv^ + 5q^v^ 
+12g2„ + lOv^q + 16w5 + y3q^) 

543 = -v^{v + lf{v + q){2q^ + 28v^ + 25v^ + 12qv + Uv'^q + llu" + 4qv^ 
+q^v^ + 2q^v + v^) 

g^i=v''{v + qf{v + lf . 
Finally, Az,s3t.j for j = 18, 19, 20 are roots of the cubic equation 

+ 53lC^ + 532^ + ff33 = 

where 

531 = -{2'iv^ + 16u^g + 16w'* + &q^v + + + q^ + 2qv^) 

932 = v'^iv + + q)iq^ + 24u^ + 26w^g + + 26v^ + 20qv^ + 5q^u^ 

+10q^v + bv^q + lOv^ + v^q^) 

g3z^~v''{v + q)\v + \f . 



The terms that enter in eq. (3.4.6) are 

^z,s3k,j — ^z,s3t,j for i = 1,2 

^z,s3k,j = Az,s3t,i+4 for 3 < j < 5 
Az,s3fcj = Az,s3ij+8 for 6 < j < 12 



(3.4.14) 



(3.4.15) 
(3.4.16) 

(3.4.17) 
(3.4.18) 

(3.4.19) 
(3.4.20) 

(3.4.21) 

(3.4.22) 
(3.4.23) 



The eight terms in the Potts model partition function for the torus strip Xz.s3t,j with 3 < j < 6 and 
10 < J < 13 do not occur for the Klein bottle strip. From our calculations of Z{G, q, v) for the Ly = 2 and 
Ly — 3 strips of the square lattice with torus and Klein bottle boundary conditions, we observe that the first 



two terms are consistent with the same generalizations (2.2.22) and (2.2.23) as for the cyclic and Mobius 
strips, although the coefficients are, in general, different. 

For the strip with torus boundary conditions the corresponding coefficients are 



cz,s3t,i =q^ - Qq^ 



- 1 



(3.4.24) 



CZ,s3t,2 



CZ,s3t,9 — ■7<-Z,s3t,3 — T^Cz,s3tA 



(3.4.25) 
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3.5 V = —1 Special Case 

For the special case v = —1, {T ~ Potts antiferromagnet) , the 20 (12) terms in the respective partition 

functions for the Ly = 3 strip of the square lattice with torus (Klein bottle) boundary conditions reduce to 
8 (5) terms entering in the chromatic polynomials for these strips. We have 

8 

P(sg,3 X m,PBCy,PBC,,q)Y,cp,s3t,j{Xp,s3t,jr (3-5.1) 

5 

P{sq,3x m,PBCy,TPBC^,q) = ^ cp,,3fej (Ap,.3fcj)" • (3.5.2) 

We order the terms Xp.sstj by decreasing degree of the associated coefficients cp,s3t,j for the torus strip, as 
polynomials in q (this differs from the order and hence labelling used in ]5l[|). These terms are 

Ap,s3t,i = -1 (3.5.3) 

AKs3t,2 = g - 1 (3.5.4) 

ARs3t,3=9-4 (3.5.5) 

Xp,s3tA = q ~2 (3.5.6) 

and 

Ap,s3t,5 = g - 5 (3.5.7) 

Ap.s3t,6 --(?'- 79+13) (3.5.8) 

\p,s3t,7 = -iq-2f (3.5.9) 

Ap..3t.8 = - 69' + 14(7 - 13 . (3.5.10) 
Thus, the correspondences for v = —1 are as follows: 

Az,s3tj ^ Ap,s3tj for j = 1,2,3 (3.5.11) 

Az,s3t,5 ~* Xp,s3t,4 (3.5.12) 

Xz,s3t,7 Xp,s3t,5 (3.5.13) 

Az,s3t,ii Ap^s3t,7 (3.5.14) 

Az,s3i,14 Xp,s3t,6 (3.5.15) 

Az,s3t,18 Xp,s3t,8 (3.5.16) 

Az,s3tj -^0 for j = 4,6,8,9, 10,12,13,15,16,17,19,20 . (3.5.17) 
The reductions for the \z.s3k,j when v = —1 follow from these and the relations ( 3.4.21 )-( t3. 4. 23| ). 
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3.6 V = Special Case 

For V = 0, i.e., infinite temperature, we have 

Xz,s3t,18 = (3.6.1) 

and 

Az,.3tj - for J + 18 (3.6.2) 
so that Z{sq, 3 x m, PBCy, {T)PBCx, g, w = 0) = g". 

4 Thermodynamic Properties of the Potts Model on Strips of the 
Square Lattice 

4.1 General 

We can use our exact results to study a number of interesting physical thermodynamic properties. Some of 
our results apply for arbitrarily great width Ly while others involve comparisons of specific widths for which 
we have carried out exact calculations. The thermodynamic properties are independent of the longitudinal 
boundary conditions (i.e., in the infinite-length direction) but do depend on the transverse boundary con- 
ditions (finite-width direction). For the strips, for any Ly no matter how large, the ferromagnet is critical 
only at T = 0, and as T — > and the correlation length ^ — > oo, the strip acts as a one-dimensional system, 
since limi^_,oo Ly/L^ — 0. In contrast, for the Potts ferromagnet on the square lattice, the phase transition 
occurs at finite temperature, at the known value Kc = ln(l + ,Jq ). Thus, studies of the thermodynamic 
behavior of the Potts model for general q ox\ LyX oo strips complement studies such as those on the approach 
to the thermodynamic limit of the Ising model on x Ly rectangular regions, in which Lx and Ly both get 
large with a fixed finite ratio Ly/L^ Q, and finite-size scahng analyses |6^, |6l| |. 

For the infinite- length limit of the Ly — Z strip with free transverse boundary conditions (i.e., the cyclic 
or Mobius strips considered here and also the free strip), the reduced free energy per site is given by 

f =\\nXz,sZcA7 (4.1.1) 

while for the infinite-length limits of the strips with periodic transverse boundary conditions (i.e. the torus 
and Klein bottle strips considered here as well as the cylindrical strip), we have, for Ly = 2, 

/= ^lnAz,s2t,5 (4.1.2) 

and, for Ly = Z, 

/= ilnAz,s3M8 • (4.1.3) 



It is straightforward to calculate from eqs. ( 4.1.1 )-( 4.1.3 ) the respective expressions for the internal energy 



U and specific heat C per site. In all cases, as noted above, the Potts ferromagnet has a T = critical 
point for general q. The Potts antiferromagnet has a. T — Q critical point if q has a certain value or values 
depending on the given strip; for the Ly = 2 strip with free or periodic transverse boundary conditions, this 
value is q = 2 p^ , while for the Ly — 3 strip with (i) free and (ii) periodic transverse boundary conditions, 
the antiferromagnet has T = critical points at (i) q — 2 and g ~ 2.34 g 13 (see eq. ( |5.1.lD below), and 
(ii) q — 2 and q — 3, respectively 
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Table 1: Some properties of the specific heat of the q-state Potts ferromagnet on strips of the square lattice with free and 
periodic transverse boundary conditions. 
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For the Ising value q — 2, since the oo hmit can be taken with even for which the cychc strip 

is bipartite, the ferromagnet and antifcrromagnet are equivalent, and hence have equivalent T — critical 



points. However, in the antiferromagnetic case, because of the noncommutativity (1.16), to avoid pathologies 
of the type discussed in detail in jl4|, one should set q = 2 first and then take the n ^ oo limit, i.e. one 
must use fnq rather than /g„. In contrast, the Ly = 3 strip with periodic transverse boundary conditions is 
not bipartite and in the q — 2 Ising case, antiferromagnetic ordering involves frustration. Nevertheless, the 
model still has a T = critical point, as is clear from the fact that the singular locus Bq runs through the 
point g = 2 in the q plane. This is somewhat analogous to the 2D Ising antifcrromagnet on the triangular 
lattice in the sense that both involve zero-temperature critical points with frustration. 

4.2 Potts Ferromagnet 

4.2.1 Free Energy, Specific Heat, and Approach to Square Lattice Limit 

In Figs. |l| and || we show plots of the specific heat for the Potts ferromagnet on the infinite- length limits 
of the Ly = 3 strips with free and periodic transverse boundary conditions. In each curve, the specific 
heat assumes a maximum value Cmax at a certain value of inverse temperature Kc,max, or equivalently, 
temperature Tc,max = J/iksKcmax)- It is of interest to inquire how these quantities depend on q and the 
transverse boundary conditions. We give this information for the results in the present paper and also our 
relevant previous studies in Table |^. In general, for the Potts ferromagnet on the infinite-length strip of 
width Ly and given transverse boundary condition, we find that, for physical q, 

Tc,max is a decreasing function of q (4-2.1) 

and 

Cmax is an increasing function of q . (4.2.2) 
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Figure 1; Specific heat (with fc^ = 1) for the Potts ferromagnet on the infinite-length, width Ly = 3 strip of the square lattice 
with free transverse boundary conditions, as a function of K = J/iksT). Going from bottom to top in order of the heights of 
the maxima, the curves are for g = 2, 3, 4. 
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Figure 2; Specific heat (with ks = 1) for the Potts ferromagnet on the infinite-length, width Ly = 3 strip of the square lattice 
with periodic transverse boundary conditions as a function oi K = J/^ksT). Going from bottom to top in order of the heights 
of the maxima, the curves are for q = 2, 3, 4. 
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The monotonicity ( [4.2. l| ) can be understood as a consequence of the fact that the specific heat, C = dU/dT, 
reaches a maximum at the temperature where there is an onset of short-range order, as reflected in the 
maximal temperature derivative of the internal energy. As q increases, there are more effective degrees of 
freedom per site, and hence it is necessary to cool the system to a lower temperature, i.e., higher value of 



K, to obtain short-range ordering. The monotonicity relation (4.2.2) shows that this onset of short-range 
order occurs more sharply as q increases. 

We may also study the dependence of Tc^max and Cmax as a function of the strip width and transverse 
boundary conditions for the Potts ferromagnet with a given value of q and investigate how Tc^max approaches 
the critical value Q 

i.e., Oc = 1 + ^/q, of the model on the infinite square lattice. As summarized in Table ^, we find that, for a 
given (physical, integral) q > 2 and for strips considered here, 

Tc.max and Cmax are increasing functions of Ly . (4.2.4) 



For the strips with free transverse boundary conditions, the monotonicity relation (4.2.4) for Tc^max can 
be explained as follows. As Ly increases, the effective coordination number, given by ( 1.28|) above, also 
increases. Since the short-range ordering is a consequence of the ferromagnetic spin-spin interactions and 
since, on average, a given spin feels the influence of ^eff of these interactions, a given degree of short-range 
order can be established at a higher temperature as Ly increases. Since, as recalled above, the maximum 
in the specific heat occurs at the temperature where this onset of short range order is taking place, the 



monotonicity relation (4.2.4) follows. Of course, in the limit Ly — > (X), the point Tc^max is a true critical 
point at which there is a nonanalyticity in the free energy, and the specific heat does not just reach a 
maximum due to short-range ordering, but diverges (logarithmically) | |59| , with the onset of long-range 
order, i.e. nonzero spontaneous magnetization. For the strips with periodic transverse boundary conditions, 
as Ly increases, there is no change in the coordination number A = 4; rather, there are actually two 
countervailing tendencies: (i) an increase in the effectiveness of short-range ordering, which, in the limit 
Ly Qo yields long-range order at sufficiently low temperature, (ii) a decrease in the finite-size effects 
involving spin-spin interactions that loop around the finite transverse periodic direction. To explain (ii), 
consider two spins which, for simplicity, are located at the same value of x but different values of y, yi and 
2/2- On an infinite square lattice, the leading contribution to the interaction of these spins would be via a 
minimum-distance path, of length r = \yi — y2\- However, for finite Ly, there is another contribution, namely 
via the path going the other way around the transverse direction, of length Ly — r. This is a finite-lattice 
artifact, the effect of which goes to zero as Ly — > oo. For the toroidal strips considered here, the effect (i) 



outweighs (ii), yielding the monotonicity relation (4.2.4). This result may be compared with the results for 



the Ising q = 2 case obtained by Ferdinand and Fisher ||58|]. These authors considered finite -size Lx X Ly 
lattices with toroidal boundary conditions and studied the limiting behavior of the specific heat as ^ cxd 
and Ly ^ oo with fixed finite nonzero ratio Ly/Lx- Among their other results, they showed that if this 
limit is taken with finite-size sections with relatively square-like aspect ratios 1/3.1 ^ Ly/L^ ^ 3.1, then 
Tc,max approaches Tc from above in the thermodynamic limit, while for more narrow rectangular shapes, 
Ly/Lx ~ 3.1 or Ly/Lx < 1/3.1, Tc^max approaches Tc from below. Although our finite-size studies of the 
q ~ 2 Potts model on toroidal strips do not, strictly speaking, fall into the category considered by Ferdinand 
and Fisher since we take Lx oo keeping Ly fixed, so that Ly/Lx is identically zero, our results may be 
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expected to correspond most closely to the extreme limit Ly/ Lx ^ in the work of [Q. Hence, one would 
expect that for q = 2 our Tc,max should approach Tc from below, and it does. It is of interest that our 
calculations generalize this monotonicity result to g > 3. 

Our calculations also show that for a q and given width Ly, the ratio Tc,max/Tc is closer to its limit 
of unity if one uses a strip with periodic rather than free transverse boundary conditions. This is easily 
explained since periodic transverse boundary conditions remove edge effects and thereby reduce finite-size 
artifacts. 

Although the (zero-field) free energy of the general g-state Potts model has not been calculated for arbi- 
trary temperature, it has been calculated in | |6^ at the phase transition temperature Tc for the fcrromagnet, 
( PI ). Of course, the free energy is also known exactly for the q — 2 Ising case. In addition to the studies 
of Ly-dependence and comparisons of Tc,max/Tc discussed in the preceding paragraphs, it is of interest to 
compare our results with a conformal field theory relation concerning finite-size effects, for the q values where 
the Potts ferromagnet has a second-order transition in 2D, namely, q — 2,3, and 4. For this purpose, we 
recall that conformal field theory methods have provided insight into the universality classes of continuous, 
second-order phase transitions and the associated critical exponents in terms of Virasoro algebras with given 



central charges and scaling dimensions [|62|-|64 . The Virasoro algebra with central extension depending on 
the central charge c is 

[L„i,L„] = (m - n)Lm+n + Y^rn{m'^ - l)S„i+n.o (4.2.5) 

For the values of q where the 2D Potts ferromagnet has a continuous transition (with infinite correlation 
length), the central charges are given by [^-|6j] (i) c = 1/2 for q = 2, (ii) c — 4/5 for q — 3, and (iii) 
c = 1 for q = 4. There is a useful relation from conformal field theory that describes the approach of the free 
energy of an infinite strip with free or periodic transverse boundary conditions to the 2D thermodynamic 
limit as oo at the critical temperature of the 2D model [p5[-p7|]: 

fstr^p,Ly = fbulk + + ^ + 0(^-3) (4.2.6) 

^y 

where fsurf. = is nonzero (zero) for free (periodic) boundary conditions in the y direction and 

r |c for PBC, 
j JLc forFBCj, 

To compare our results with this relation, we use as input the critical values of the Potts model free energy 
on the square lattice. For q ~ 2, from the Onsager solution one has 

fpotts,q=2,crtt. = J In 2 + ln(l + V2) + -Gc^ 1.811068985... (4.2.8) 

2 TT 



where Gc is the Catalan constant 



For q = 3 and g = 4, one has 



= XI T^TT-^ = 0-91596559... (4.2.9) 

k=0 



(2fc+l)2 



11 r- iGr 

/potta,q=3,c„t. =ln2+-ln3 + -ln(2 + V3) + —^ = 2.07018716... (4.2.10) 
2 6 on 



and 

TPoUs, q—4,crit. 



fpoUs,,=4,cr^t. = -31n2 + 41n( ^ 2.25952475... (4.2.11) 
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Table 2: Comparison of exact partition function calculations with conformal field theory finite-size relation for the infinite- 
length, width Ly = 3 strip with periodic transverse boundary conditions. See text for further details. 



q 


RHS 


f strip (-^c) 


f strip i^^C .max ) 


2 


1.840 


1.8425 


1.91 


3 


2.117 


2.121 


2.18 


4 


2.318 


2.324 


2.37 



where T{z) is the Euler Gamma function. As an illustration, we consider the width Ly = 3 strip with 
periodic transverse boundary conditions. In Table |^ we list the values of the right-hand side of (4.2.6) and 
the evaluation of fstrip at the value T — where the model is critical on the 2D lattice. Of course, for 
any finite value of Ly, fstrip is not critical at this value of temperature and, for comparison, we also list 
the comparison with fstrip evaluated at T = Tc,max- Clearly, as Ly oo, Tc,max Tc, so that these are 



equivalent in this limit. One sees that eq. ( 4.2.6 ) is reasonably well satisfied, even for the modest value of 



the width Ly ^ 3. 



4.2.2 T ^ Limit 



We next derive some general results concerning the low-temperature behavior of the (reduced) free energy 
/ and, from this, the internal energy U and specific heat C, for the g-state Potts ferromagnet on an infinite- 
length strip of the square lattice with arbitrary longitudinal boundary conditions and specified transverse 
boundary conditions. For technical convenience, we use free longitudinal boundary conditions for our calcu- 
lations, assume a large length L^; when calculating the partition function, and then, as usual, take L^ ^ oo 
when calculating the free energy (per site). At T — 0, Z — qa'^'^'^\ where, as defined above, e(G) denotes 
the number of edges (bonds) on the strip graph. Now consider finite-temperature corrections to this result. 
If the lattice were of dimensionality ci > 2, the leading order term always arises from a configuration in 
which one changes one spin from its preferred value to one of the (g — 1) other values, and so forth for 
successive corrections. However, the leading corrections do not, in general, arise in this simple a manner for 
the finite- width strips, especially in the case where the transverse boundary conditions are free. Let us begin 
with these. 

For Ly = \, the leading low-temperature correction arises from changing all of the spins to the right, say, 
of a given point from their previous value to one of the (q — 1) other values. Since this can be done at any 
of Lx locations, we obtain f = K + {q — l)e^^ + where ... here and below indicate subdominant terms 
in the T — > limit, and hence U = J[— 1 + (g — l)e~^] and C — ksK^iq — l)e^^ . Of course, in this Ly = 1 
case, the exact expressions for /, U , and C are simple: / = ln(e'^ + q — \), U = —Je^/{e^ + q — 1), and 
C = ksK^q - l)e^/{e'^ + q- l)\ 

For Ly — 2, the leading correction to the low-temperature limit for Z is obtained by changing all of 
the spins to the right of a given transverse slice from their original value to one among the {q — 1) other 
possibilities. This yields the expansion 

fFBC,,L,=2 = lK + ^{q- l)e-2^^ + ... (4.2.12) 



and hence 



-^ + (g-l)e-2^ + ... (4.2.13) 
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and 



-2K 



CFBC,.L,=2 = 2kBK\q-l)e- 

I, in agreement with a direct derivation from the exact results obtained therein. 



(4.2.14) 



as given m 

For Ly = 3, there are two different types of leading corrections to the low-temperature limit for 
obtained by (i) changing all of the spins to the right of a given transverse slice from their original value to 
one among the {q — 1) other possibilities, and (ii) changing a spin on the upper or lower boundary from its 
original value to one among the [q — 1) others. These both involve a suppression factor e~^^ , and one gets 
the expansion 



fpBCy 



-K+{q-l)e- 



-3K 



and hence 



and 



FBCy,L. 



-'- + 3{q-l)e-'^< 



-3K 



(4.2.15) 
(4.2.16) 
(4.2.17) 



CFBCy..L,=3^9kBK'{q-l)e- 

We have that this is consistent with a direct derivation from the exact result ([4.1.l| ). 

For greater widths, Ly > 4, the leading low-temperature correction arises from changing one of the spins 
on the upper or lower horizontal edges from their original values to one of the other (q — 1) possible values. 
This makes the dominant contribution because of the fact that these edge spins have a lower coordination 
number (degree), namely A = 3, than the spins in the interior of the strip. We have 



FBCy.Ly>4 



^eff 



K + —{q~l)e-'^'^ + 



where ^eff was given in eq. ( |l.2§| ) above, and hence 



UpBCy 



,L„>4 



L. 



l)e 



-3A' 



and 
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CFBCy,Ly>i = —ksK^iq - l)e 



-3K 



(4.2.18) 



(4.2.19) 



(4.2.20) 



The essential zeros in U and C are typical for a spin model at its lower critical dimensionality and mean that 
the low-temperature series expansion has zero radius of convergence. Physically, this reflects the qualitatively 
different behavior at T = 0, where there is saturated magnetization, and at any nonzero temperature, 
regardless of how small, where this magnetization vanishes identically. 

Now consider a strip of the square lattice with free transverse boundary conditions and arbitrarily great 
width Ly > 3. Let the free energy be denoted as 



1 



f ^ -j— In Al^^pm 



(4.2.21) 



where \Ly,PM is the term in (1.9) that is dominant for physical temperature (paramagnetic phase). Then 



for Ly ^ 3 this has the low-temperature expansion 



Ar,= 



3,PM 



„5K 



l + 3(g- 1). 



-3K 



(4.2.22) 
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i.e., in the notation of eq. ( |1.19| ), 

Ai -3,PM,u = 1 



3(g-l) 



for 



and for Ly > A, 



^ p(Ae///2)L^/f 



^Ly,PM 

and hence, in the notation of eq. ( 1.19| ), 



1 + 2(9- l)e- 



-3K 



for 



. 



(4.2.23) 
(4.2.24) 

(4.2.25) 



Since the singular locus B is defined in the neighborhood of the point u = by the degeneracy of magnitudes 
|ALy,PAf| — \^Ly,j\, where ^Ly,j is another leading term at this point, it follows that this locus does not 
involve an arbitrarily large number of curves crossing in an intersection point at u = 0, and the corresponding 
complex-temperature partition function zeros (Fisher zeros ) do not become dense in the vicinity of the 
origin as Ly —^ 00. These results supercede a conjecture given in 

We next derive corresponding low-temperature expansions for strips with periodic transverse boundary 
conditions. In the cases Ly = 2 and Ly ~ 3, the leading correction term arises from changing all of the spins 
to the right of a transverse slice from their original values to one of {q — 1) other values. This yields 



fPBCy,Ly=2K+^iq-l)e 



-2 + {q-l)e 



-LyK 



CpBCy,Ly = LyksK^iq - l)e 



-LyK 



(4.2.26) 

(4.2.27) 
(4.2.28) 



For Ly = 4, these spin configurations involving a seam make a contribution equal to those from configurations 
in which one changes one of the spins on a transverse slice to one of [q — 1) other values, so that 



fpBCy.Ly=i = 2K+-{q-l)e 



-4K 



UpBCy 



J 



(4.2.29) 

(4.2.30) 
(4.2.31) 



-2 + 5(g-l)e-4^ + . 

CpBCy,Ly=4 - 20kBK\q - l)e-4^ + . 

Finally, for Ly > 5, the leading correction term in the low-temperature limit arises in the same way as on an 
infinite 2D square lattice with toroidal boundary conditions, namely via the change of a single spin to one of 
{q — 1) other values. This yields expansions whose leading term is the same as on the infinite square lattice: 

fpBCy.Ly>5 = 2K+{q-l)e-^'' + ... 



lPBCy,Ly>5 
UpBCy,Ly>5 — J 



-2 + 4(<z-l)e-*^ + 



CpBCy,Ly>5 - l&kBK^{q-l)e 



+ .. 



(4.2.32) 
(4.2.33) 
(4.2.34) 



Next consider a strip of the square lattice with periodic transverse boundary conditions and arbitrarily 
great width Ly > 5. Let the free energy be denoted as in (4.2.21). Then we have the low-temperature 
expansion 



^Ly,PM — e " 
i.e., in the notation of eq. ( |l.l9| ), 

^Ly,PM,U 



l + Ly{q-l)e 



-4A' 



l + Lyiq-l)u^ + 



for u — !■ 



(4.2.35) 
(4.2.36) 
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Table 3: Some properties of the specific fieat of the q-state Potts antiferromagnet on strips of the square lattice with free and 
periodic transverse boundary conditions. Since for q > 3 this model has no finite-temperature phase transition on the square 
lattice, the last column is left blank. 



q 


BCy 


Ly 


,max 


Cmax / 




2 


FBCy 


1 


2.40 


0.439 


0.367 


2 


FBCy 


2 


1.33 


0.631 


0.663 


2 


FBCy 


3 


1.17 


0.755 


0.753 


2 


PBCy 


2 


1.00 


0.529 


0.8805 


2 


PBCy 


3 


1.88 


0.366 


0.469 


3 


FBCy 


1 


2.23 


0.241 




3 


FBCy 


2 


1.785 


0.3325 




3 


FBCy 


3 


1.70 


0.369 




3 


PBCy 


2 


1.33 


0.287 




3 


PBCy 


3 


2.62 


0.498 




4 


FBCy 


1 


2.16 


0.166 




4 


FBCy 


2 


1.96 


0.238 




4 


FBCy 


3 


1.92 


0.263 




4 


PBCy 


2 


1.49 


0.208 




4 


PBCy 


3 


2.16 


0.346 





4.3 Potts Antiferromagnet 

4.3.1 Free Energy, Specific Heat, and Approacli to Square Lattice Limit 

In the case of the Potts antiferromagnet, as discussed in [|l^ |l^, there can be nonanalyticities at sufficiently 
small non-integral q and temperature, but these do not represent physical phase transitions and involve a 
number of pathologies such as non-existence of a thermodynamic limit independent of longitudinal boundary 
conditions, negative specific heats, and lack of a Gibbs measure. For the physical cases, the analyticity of 
the free energy is a consequence of the general theorem that one-dimensional and quasi-one-dimensional spin 
systems with short-range spin-spin interactions do not have a finite-temperature phase transition, which, in 
turn, is proved by an elementary application of a Peierls argument. This analyticity property is equivalent 
to the property that the singular locus Bu does not cross the positive real u axis. 

In Figs. ^ and ^ we show plots of the specific heat for the Potts antiferromagnet on the infinite-length 
limits of the Ly — i strips with free and periodic transverse boundary conditions. Using our exact results, 
we have calculated similar plots for Ly = 2, but to save space we do not show them here. In Table ^ we list 
some comparative information on Kc,max and Cmax from our previous and current exact calculations. In 
contrast to the situation with the ferromagnet, l-fiTcmaxI does not have a uniform behavior with respect to 5; 
for example, for the strips with free transverse boundary conditions, for Lj, = 1, it decreases with increasing 
q while for Ly ~ 2 and Ly — i it increases with increasing q. For the strips with periodic transverse boundary 
conditions, lii'c.maa;! increases with increasing q for Ly = 2 but is not a monotonic function of q for Ly = 3. 
Note that for the q = 2 Ising value on the Ly — 3 strip with periodic boundary conditions the frustration 
affects the value of Kc,max', this would also be true on wider strips of this type with odd Ly. For fixed q, 
as a function of Ly, we find that for the strips with free transverse boundary conditions, \Kc,max\ decreases 
and Cmax increases as Ly increases; for the strips with periodic boundary conditions these trends also apply 
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Figure 3: Specific heat (with ks = 1) for tlic Potts antiferromagnet on the infinite-length, width Ly = 3 strip of the square 
lattice with free transverse boundary conditions, as a function oi K = J/iksT). Going from top to bottom in order of the 
heights of the maxima, the curves are for g = 2, 3, 4. 
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Figure 4: Specific heat (with ks = 1) for the Potts antiferromagnet on the infinite-length, width Ly = 3 strip of the square 
lattice, as a function of K = J/iksT). Going from top to bottom in order of the heights of the maxima, the curves are for 
9 = 3,2,4. 
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for g > 3, but for q = 2, one encounters the effects of frustration for Lj^ = 3 and, more generally, for odd 
Ly. Since the g = 3 Potts antiferromagnet has a zero-temperature critical point on the square lattice, one 



can make use of the relation (4.2.6) from conformal field theory. In this case the model exhibits ground 



state entropy without frustration, and, since Yitcit^o PU — 0, the general relation / = ~[3U + S/ks (i.e. 



F = U — TS) reduces to / = Sa/ks- Hence, eq. (4.2.6) becomes 



Ssurf. A _3 

Sstrip,Ly = Sbulk H J — - + + C>{Ly ) . (4.3.4) 

Using the exact calculation of the ground state entropy for the Ly = A strip with cylindrical boundary 
conditions |Q, it was noted earlier that this relation is satisfied quite well for the known value of the 
central charge, c — 1. 

4.3.2 T ^ Limit 

As regards the low-temperature expansions of the internal energy and specific heat for the antiferromagnetic 
Potts model on these strips, these depend on g in a more complicated manner than the analogous expansions 
for the ferromagnetic case, as was already clear from eqs. (6.35)-(6.37) and (6.39) in iQ. For the Ising value 
q = 2 on strips of the square lattice with free transverse boundary conditions, the limit ^ oo can always 
be taken such that the graphs are bipartite, and hence the correction terms in U and C are obtained from 
those for the ferromagnet by the replacement K — > —K. However, for other values of g, the dependence is, 
in general, different. For example, for Ly — 2, C has an exponential factor e^^ if g = 2 but e^' if g 7^ 2 [^ . 

5 Locus B for the Ly = 3 Strip with CycUc or Mobius Boundary 
Conditions 

5.1 B in the q Plane 

The singular locus B in the q plane was given in for the special case a = 0, i.e., the zero-temperature 
Potts antiferromagnet. This locus separates the q plane into seven regions. The first of these is Ri, which 
includes the real intervals q < and q > qc, where 

qc ^ 2.33654.. for a = 0, G = sq, Ly = 3, FBCy, {T)PBC^ . (5.4.4) 

The other regions are i?2, including the real interval 2 < q < qc', R3 including the real interval < g < 2; 
two complex-conjugate (c.c.) phases i?4,i?| centered at g ~ 2.4 ± 0.9i; and two additional quite narrow, 
sliver-like c.c. phases i?5, centered at q ~ 4.93 ± 4.70i. 4n region i?i the dominant A is, with our current 
ordering of terms, Xp^s3c.9j so that for either cyclic or Mobius longitudinal boundary conditions, 

W{sq, 3 X 00, FBCy, {T)PBC,, q) = (Ap,,3c,9)'/' for q e Ri . (5.4.2) 
In region i?2, the dominant A is q — 4, so 

\W{sq,Sx 00, FBCy, {T)PBC^,q)\ = \q- for q £ R2 ■ (5.4.3) 
Recall that for regions other than Ri , only the magnitude of W can be determined unambiguously [53 . In 



region i?3, the dominant A is the maximal root of the cubic ( 2.3.8 ). In the complex conjugate pairs of regions 
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Figure 5: Singular locus B in the q plane for the free energy of the Potts antifcrromagnct for the temperature given by a = 0.1, 
i.e., K = — In 10, on the width Ly = 3, infinite-length strip of the square lattice with cyclic or Mobius boundary conditions. 
Partition function zeros for the cyclic strip with length Lx = 20 and thus n = 60 vertices are shown for comparison. 
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Figure 6: Singular locus B in the q plane for the free energy of the Potts antifcrromagnct for the temperature given by a = 0.5, 
i.e., i<r = — ln2, on the width Ly = 3, infinite-length strip of the square lattice with cyclic or Mobius boundary conditions. 
Partition function zeros for the cyclic strip with length Lx = 20 and thus n = 60 vertices are shown for comparison. 
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(i?4, -R4) and (i?5, R^) the dominant A's are the other two roots of the cubic ( 2.3.8 ). In each of these regions, 
\W\ = |Ap.s3cj-,(iom. r^^, where \p^s3cj,dom. IS the dominant root in the respective region. 

In Figs. 1^ and ^ we show plots of the locus Bq for the Potts antiferromagnet at the values of temperature 
given by a = 0.1, i.e., _ftr = — In 10, and a = 0.5, i.e., K ^ — In 2. For comparison, typical sets of partition 
function zeros in the q plane arc shown for long finite strips. As a increases from to 1, i.e., the temperature 
increases from to infinity for the Potts antiferromagnet, the boundary Bq contracts and finally shrinks to 
a point at the origin in the q plane. In the interval < a < acre, where acre = 0.181720.., the rightmost 
part of the boundary and hence also qc sweep to the left, while the crossing at q = 2 remains; at a = ac e, 
the rightmost part of the boundary coincides with q ~2, i.e., qc — 2. As a increases further in the interval 
acre a < 1, there are only two regions, Ri and R^, that include real intervals. For this full range of 
temperatures in the Potts antiferromagnet, < a < 1, the boundary B always passes through q = 0. 



Im(q) - 




-3 I ' ' ' ' ' ' ' ' ' ' ' 1 

-5 -4 -3 -2-10 1 

Re(q) 



Fi gure 7'. Singular locus 3 in the q plane for the free energy of the Potts ferromagnet for the temperature given by a = 2, i.e., 
A" = In2, on the width Ly = 3, infinite-length strip of the square lattice with cyclic or Mobius boundary conditions. Partition 
function zeros for the cyclic strip with length Lx = 20 and thus n = 60 vertices are shown for comparison. 

In Fig. we show the locus Bq for a typical ferromagnetic value, a = 2, i.e. if = In 2. One can discern 
very small prongs extending in northwest and southwest directions from the main curve at g ~ — 2.74± 1.88i, 
and a short line segment on the negative real q axis. Two general features are that the locus B (i) passes 
through q — Q and (ii) does not cross the positive real. 
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5.2 B in the u Plane 



lm(u) 



-2 -1 



1 

Re(u) 



Figure 8: Singular locus B in the u plane for the free energy of the Potts model for g = 2, on the width Ly = 3, infinite-length 
strip of the square lattice with cyclic or Mobius boundary conditions. Partition function zeros for the cyclic strip with length 
Lx = 20 and hence e = 100 are shown for comparison. 

In Figs. I^pl] we show plots of the singular locus B for the Potts model free energy for q — 2,3, and 10, 
respectively, in the infinite-length limit of the Ly — 3 strip graph with cyclic or Mobius boundary conditions. 
Zeros of the partition function are shown for comparison. There are several interesting features of these 
plots. In each, there are six curves forming three branches on B that intersect at the origin, u = 0, at 
the angles ±7r/6, ±7r/2, and ±57r/6. In general, the density of complex-temperature Fisher zeros along the 
curves comprising B in the vicinity of a generic singular point Us behaves as |7l|, |7^ 

g^\u-us\^-"^ (5.2.1) 

where as (a^) denotes the corresponding specific heat exponent for the approach to Us from within the CTE 
PM (FM) phase. Thus, for a continuous, second-order transition, with < 1, this density vanishes as one 
approaches the critical point along B. The specific heat for these infinite-length, finite-width strips has 
an essential zero at T = which, if expressed in terms of an algebraic specific exponent a, corresponds to 
= — ooatu = Us = 0. Substituting this into ( ^.2.1 ), it follows that the density vanishes rapidly as one 



a 



approaches the origin u = along the curves forming Bu- This is evident in Figs. PjTl|. As we have discussed 
in our earlier works on B in the q and u plane, any mapping of the regions in these respective planes is done 



34 



3 



2 



1 



lm(u) 



-1 



-2 



-3 




-2-10 1 
Re(u) 

Figure 9: Singular locus B in the u plane for the free energy of the Potts model for q = 3, on the width Ly = 3, infinite-length 
strip of the square lattice with cyclic or Mobius boundary conditions. Partition function zeros for the cyclic strip with length 
Lx = 20 and hence e = 100 are shown for compaxison. 
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Figure 10: Singular locus B in the a = u^^ plane for the free energy of the Potts model for q = 3, on the width Ly = 3, 
infinite-length strip of the square lattice with cyclic or Mobius boundary conditions. Partition function zeros for the cyclic strip 
with length Lx = 20 and hence e = 100 are shown for comparison. 
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Figure 11: Singular locus B in the u plane for the free energy of the Potts model for q = 10, on the width Ly = 3, infinite- 
length strip of the square lattice with cyclic or Mobius boundary conditions. Partition function zeros for the cyclic strip with 
length Lx = 20 and hence e = 100 axe shown for comparison. 
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with finite resolution and does not exclude the possibility of extremely small regions. That such regions can 
occur was evident, e.g., in Q. 

In earlier work ||7^, 17 1, it was shown that although the physical thermodynamic properties of 

a discrete spin model are, in general, different in quasi-lD systems such as infinite- width, finite- width strips, 
and in higher dimensions, nevertheless, exact solutions for Bu in quasi-lD systems can give insight into Bu in 
2D. This was shown, in particular, for the q = 2 Ising special case of the Potts model, where the comparison 
can be made rigorously since the model is exactly solvable in 2D. Our current results extend this comparison 
to greater width. For example, in the complex-temperature phase diagram of the Ising model on the square 
lattice, the singular locus B forms two circles u = ±1 + v^e**^, < (/) < 27r This locus thus exhibits 

complex-conjugate multiple points at u = ±i, where four curves forming two branches on B intersect. This 
was also true of the Ising model on the infinite-length, width Ly = 2 cyclic or Mobius strips of the square 
lattice 

For the q — 2 (Ising) case shown in Fig. ^, the locus B is invariant under the inversion mapping 
u l/u = a. This is a consequence of the fact that the infinite-length limit can be taken on a sequence of 
bipartite graphs. It follows that the six curves extending outward to infinity in the u plane cross at the origin 
of the a plane. Thus, the points u = and a — are multiple points on B. There are also complex-conjugate 
multiple points at u = ±i, where again six curves forming three branches of B cross each other. Finally, 
there is another multiple point at u = — 1 where four curves forming two branches intersect. Although the 
locus B does not appear to contain any arc endpoints for q — 2, such endpoints are present for the other 
illustrative values of q, namely q — i and q = 10. 

In the q — 3 case, the intersection points of the various curves include the complex-conjugate pair 
u ~ q±2tti/3 .^gjj g negative real point, u = —1/2. The intersection points include both those where 
several curves cross and those where curves come together in a T-type intersection point. A previous 
analytic discussion of such T-type intersection points in complex-temperature phase diagrams was given in 
|f73| . Besides the origin, u = 0, the locus B crosses the real u axis also at the three points u — —2, u ~ —1.37, 
and u ~ —1.21. We recall that for g = 3 Potts model on the infinite-length, width Ly = 2 cyclic or Mobius 
strip of the square lattice, B also included multiple points at u = e^^'^^^^ and crossed the real axis at m = —2 
and u = —1/2 (as well as the origin) ||l^. These results for the Ly = 2 and Ly = 3 cyclic/Mobius strips 
suggest that it is likely that for the infinite 2D square lattice, B passes through the points u — e^'^'^^^^. 
Calculations of Fisher zeros on finite lattices show considerable scatter in the region Re{a) < |7^, |79| 
but are consistent with this suggestion. The crossings at u = —1/2, i.e., a ~ —2, for B on these Ly = 2,3 
strips are of interest since, given that the g = 3 Potts antiferromagnet has a zero-temperature critical point 
on the square lattice, so that the singular locus B passes through a = 0, it follows by duality that this locus 
also passes through a = —2 fl^ . Since some features of B are clearer when it is shown in the a = u^^ plane, 
we do this in Fig. |l^. 

The plot of 6 for g = 10 plot in Fig. shows another interesting feature: a portion of this locus, and 
the associated partition function zeros exhibit an approximately circular shape. Of course the actual locus 
is considerably more complicated than a circle. However, in this context, one may recall that in the q oo 
limit, the complex-temperature phase boundary is |C| = 1, where C = (a — 1)/%/? (T^j which, for large q, 
yields the locus |u| 1/v^ in the u plane. Thus, aside from the features of B that reflect the quasi-lD 
nature of the strip, such as the curves passing through u — 0, the width = 3 is large enough so that one 
can begin to see this feature of the g — s- cx) limit. 
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6 Locus B for Ly = 2 Strip with Torus/Klein Bottle Boundary 
Conditions 

We find that the locus Bq for the (infinite- length limit of the) Ly = 2 strip with torus or Klein bottle boundary 
conditions is generally similar to that for the same strip with cyclic or Mobius boundary conditions [p^ . For 
the finite-temperature ferromagnet this is again true. In particular, Bq for the ferromagnet does not cross 
the positive real q axis. 

In Figs. ^ and |l^ we plot Bu for g = 2 and q = 3. These loci contain four curves forming two 
branches that cross at the origin, with angles ±7r/4 and ±37r/4. Comparing the locus for g = 2 plot with 
the corresponding locus in the case of the infinite-length Ly = 2 cyclic or Mobius strip of the square lattice 
|l5|| we observe that here the locus is symmetric under the map u — > — u, whereas for the cyclic or Mobius 
strip it is not. The reason for this is that the torus or Klein bottle strip has an even coordination number, 
A = 4, whereas the cyclic or Mobius strip has an odd coordination number, A = 3 (both are A-regular 
graphs). The g = 3 locus contains intersection points at u — e±27rj/3 ^^^1 as two pairs of T (equivalently, 
Y) type intersection points, and crosses the negative real axis at u = —1 and u — —2. 



Im(u) 



-2 



-3 -2 -1 
Re(u) 



Figure 12: Singular locus B in the u plane for the free energy of the Potts model for g = 2, on the width Ly = 2, infinite-length 
strip of the square lattice with torus or Klein bottle boundary conditions. Partition function zeros for the toroidal strip with 
length Lx = 20 and hence e = 80 are shown for comparison. 



39 



Im(u) 



-1 



-2 



-1 
Re(u) 



Figure 13: Singular locus B in the u plane for the free energy of the Potts model for q = 3, on the width Ly = 2, infinite-length 
strip of the square lattice with torus or Klein bottle boundary conditions. Partition function zeros for the toroidal strip with 
length Lx = 20 (e = 80) are shown for comparison. 



40 



7 Locus B for Ly = 3 Strip with Torus/Klein Bottle Boundary 
Conditions 

7.1 B in the q Plane 



2 



lm(q) - 



-1 



-2 



Re(q) 

Figure 14: Singular locus B in the q plane for the free energy of the zero-temperature Potts antiferromagnet on the Ly = 3 
strip of the square lattice with torus or Klein bottle boundary conditions. Chromatic zeros for the toroidal strip with length 
Lx = 30 and thus n = 90 vertices are shown for comparison. 

For comparison with our finite-temperature results, we show in Fig. |lj the locus B for the infinite-length 
limit of the Ly = 3 strip with torus or Klein bottle boundary conditions in the q plane for a — 0, i.e., 
the zero-temperature Potts antiferromagnet, on the Ly = 3 strip of the square lattice with torus or Klein 
bottle boundary conditions This locus separates the q plane into three regions: Ri including the real 
intervals q > qc, where qc — 3, and q < and extending to arbitrarily large \q\; i?2, including the real interval 
2 < q < 3, and (iii) R3, including the real interval < g < 2. Thus, Bq crosses the real q axis at q — 0, 2, 
and 3. For the same model at the finite values of temperature given by a = 0.1 and a — 0.5, we show Bq 
in Figs. |l^ and |l^. As a increases from to 1, i.e., the temperature increases from to infinity for the 
Potts antiferromagnet, the boundary Bq contracts to a point at the origin in the q plane. In the interval 
< a < Ocr, where Ucr = 0.3488338..., which is a root of the equation 



2z 
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13^' 



2lz 
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Figure 15: Singular locus B in the q plane for the free energy of the Potts antiferromagnet for the temperature given by a = 0.1, 
i.e., K = — In 10, on the width Ly = 3, infinite-length strip of the square lattice with torus or Klein bottle boundary conditions. 
Partition function zeros for the toroidal strip with length Lx = 20 and thus n = 60 vertices are shown for compajrison. 
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Figure 16: Singular locus B in the q plane for the free energy of the Potts antiferromagnet for the temperature given by a = 0.5, 
i.e., = — ln2, on the width Ly = 3, infinite-length strip of the square lattice with torus or Klein bottle boundary conditions. 
Partition function zeros for the toroidal strip with length Lx = 20 and thus n = 60 vertices are shown for compajrison. 
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-llz^ + 44z^ - 70z^ + 35z2 - 17z + 4 = 



(7.1.1) 



The boundary between regions i?2 and i?3 remains fixed at q = 2 and the right-most part of tlie boundary, 
passing through qc and separating Ri and i?2 moves leftward, and crosses q = 2 at a = Ocr. In the interval 
acr < a < 1, there are only two regions, Ri and i^s, that contain intervals of the real axis, and continues 
to decrease from 2 to 0. 

In Fig. ^ we show the locus Bq for a typical ferromagnetic value, a = t; + 1 = 2, i.e. if = In 2. As 
is evident in this figure, B includes a complex-conjugate pair of prongs extending outward from the main 
curve. 
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Figure 17: Singular locus B in the q plane for the free energy of the Potts ferromagnet for the temperature given by a = 2, 
i.e., K = \n2, on the width Ly = 3, infinite-length strip of the square lattice with torus or Klein bottle boundary conditions. 
Partition function zeros for the toroidal strip with length = 20 and thus n = 60 vertices are shown for comparison. 

7.2 B in the u Plane 

We show here plots of the singular locus B and partition function zeros in the u = 1/a plane for several 
illustrative values of q, namely, q = 2, q = 3, and q = 10 in Figs, [l^ pO| . 

Some of the features of these plots are similar to those discussed above for the infinite-length Ly = 3 strips 
of the cyclic or Mobius strips, such as the sixfold crossing of curves at u = with density of Fisher zeros 
rapidly vanishing near this point, and the additional multiple points at u = ±i for q = 2 and at u — e±27ri/3 
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Figure 18: Singular locus B in the u plane for the free energy of the Potts model for g = 2, on the width Ly = 3, iiifiiiite-leiigth 



strip of the square lattice with torus or Klein bottle boundary conditions, 
length Lx = 15 and hence e = 90 are shown for compajrison. 



Partition function zeros for the toroidal strip with 
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Figure 19: Singular locus B in the u plane for the free energy of the Potts model for q = 3, on the width Ly = 3, iiifiiiite-leiigth 
strip of the square lattice with torus or Klein bottle boundary conditions. Partition function zeros for the toroidal strip with 
length Lx = 20 and hence e = 120 are shown for comparison. 
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Figure 20: Singular locus B in the u plane for the free energy of the Potts model for q = 10, on the width Ly = 3, infinite- 
length strip of the square lattice with torus or Klein bottle boundary conditions. Partition function zeros for the toroidal strip 
with length Lx = 20 and hence e = 120 are shown for compajrison. 
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for g = 3. Note that in the present case, for q = 3, B crosses the real u axis at the points u = — 1 and 
u = —(1 + •\/3/)2 ~ —1.366..., in addition to the origin. In contrast to the Ly = 3 cychc or Mobius strips, 
where B is compact in the u plane, here, for the torus or Klein bottle strips, two complex-conjugate curves on 
B extend to complex infinity in the u plane, i.e., pass through the origin of the a plane. This is equivalent to 
the property shown in |Q that B for the zero-temperature (i.e., a = 0) Potts antiferromagnet passes through 
q — Qc — 3. As was noted in |5ll , this is interesting, since it indicates that the q = 3 Potts antiferromagnet 
has a zero-temperature critical point on the infinite-length, width Ly = 3 strips of the square lattice with 
torus or Klein bottle boundary conditions. As we have discussed before, a convenient feature of strips with 
periodic longitudinal boundary conditions is that for the g-state Potts antiferromagnet with a value of q = qc 
where the model has a zero-temperature critical point, this is signalled explicitly by the properties that Ba 
passes through the origin of the a plane and Bq passes through q = qc for a — 0. In contrast, if one uses 
strips with free longitudinal boundary conditions, this connection is not evident in general. Thus, if one uses 
a cylindrical rather than toroidal or Klein bottle strip with Ly = 3, then the boundary Bg does not pass 
through q = 3 H. 

8 Conclusions 

In this paper we have presented and discussed exact calculations of the partition function of the g-state Potts 
model for general q and temperature on strips of the square lattice of width Ly = 3 vertices and arbitrary 
length Lx with longitudinal boundary conditions of the cyclic, Mobius, torus, and Klein bottle types. In 
the infinite-length limit we analyzed the resultant thermodynamic properties and also derived a number of 
low-temperature expansions, including results valid for arbitrarily wide strips with free or periodic transverse 
boundary conditions. A number of interesting results were given for the continuous locus B where the free 
energy is singular in both the q plane and temperature variable plane. We also discussed the related Tutte 
polynomials. 
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9 Appendix 

9.1 Equations for Xz,s3, 



The terms Xz,s3c,j for 10 < j < 15 in eq. (2.2.1) and (2.2.2) are roots of the sixth degree equation 



where 



e + hl^" + h2C + /63^' + + /65e + ./66 = 

/ei = -w(7f^ + 18v^ + v^q + IQqv + 2q^ + v"^) 
/g2 = v^{q^ + Vlq^v + 2q^v^ + q + 14u^g + 2v^(^ + 2w^ + 19^;^ 
-|-78u^ + 64t;^ -I- 109gu^ -I- h^c^v^ -I- 66w^g -t- 201;^^^) 



(9.1.1) 
(9.1.2) 

(9.1.3) 



/63 = -v^(y ^ q){2q^ + 21q^v + I3q\^ + 8v^q + 62v'^q + 19v'^q^ + 16w^ 
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+79v^ + I26v^ + + q^v + q^v^ + 2v^q^ + 168v^ + 168qv^ + 9lq^v^ 

+l^2v^ q + nv^q^) (9.1.4) 

/64 = v^{v + l){v + qf{q^ + liq^v + ^q^v^ + 2v^q + 2Qv^q + Tu^'q^ + Zv^ 

+27v^ + 78v^ + q^v^ + 77v^ + llOqv^ + 56q^v^ + SOw^g + 33v^q^) (9.1.5) 

/eg = -w9(w + g)4(„ + 1)2(18^2 + llqv + 7v^q + + qv^ + Uv^ + q^v + 2q^) (9.1.6) 

h^^v'\v + qf{v + lf. (9.1.7) 
The terms \z,s3c,j for 17 < j < 20 are roots of the quartic equation 

+ hi^ + /42e' + AaC + /44 = (9.1.8) 

where 

/4i = -(15w^ + Gu"* + qv^ + 12v2q + 5^2^ + + w^) (9.1.9) 
/42 = w^(w + q){2q^ + 32u3 + SQw^g + + 30u'' + 8q\2 + 13^2^ ^ 4^4^ 

+ g^w + 21qw3 + lOu^) (9.1.10) 
/43 = -w4(„ + l)(t; + qf[q^ + 15„3 4- i9„2q ^ow'* + _,. 9^2^ 

+2A + v^g^ + 10qw^ + 2u5) (9.1.11) 

9.2 Relations Between Potts Partition Function and Tutte Polynomial 

The formulas relating the Potts model partition function Z{G^ q, v) and the Tutte polynomial T{G, x, y) were 
given in jl4j and hence we shall be brief here. The Tutte polynomial of G, T{G, x, y), is given by (7j-|^ 

T{G, x,y)=J2{x- l)^(G')-MG)(y _ i)c(G') (g 2.1) 

G'CG 

where k{G'), e(G"), and n{G') — n{G) denote the number of components, edges, and vertices of G', and 

c(G') = e(G') + fc(G') - n(G') (9.2.2) 
is the number of independent circuits in G'. For the graphs of interest here, k{G) = 1. Now let 

x=l + - (9.2.3) 

V 

and 

y^a^v + l (9.2.4) 
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so that 



q={x-l)iy-l) 



Then 



Z{G,q,v) - (x- l)^(«)(y- ir(«)T(G,.T,2;) 
For a planar graph G the Tutte polynomial satisfies the duality relation 



(9.2.5) 
(9.2.6) 

(9.2.7) 



where G* is the (planar) dual to G. As discussed in [Q , the Tutte polynomial for recursively defined graphs 
comprised of m repetitions of some subgraph has the form 



T{G„i,x,y)= ^ ct,Gj(At,Gj)^ 



(9.2.8) 



9.3 Cyclic and Mobius Ly = 3 Strips of the Square Lattice 

We calculate 

20 

T(s9,3 X m,FBCy,PBC^,x,y) = J] ct,.3c,, (At,s3cj)" (9.3.1) 

(do not confuse the a: and y referring to the longitudinal and transverse directions with the x and y variables 
defined in (|9.2.3|) and (|9.2.4|)) and 



where 



T{sq, 3 X m, 

FBCy,TPBC,,x,y) = J2 



At,s3c,1 = 1 
At,s3c,2 = X 



A 



r,s3c,(3,4) 



a; + y + 3 ± a/ (a; + y + 3)2 - ixy 



and 



1 



At,s3c,(5,6) = ^ x + y + l± v/(a; + y + 1)2 - 4xy . 
Among the AT,s3c,j's that depend on x and y, four are symmetric functions of these variables: 

>^T.s3c.jix,y) = XT.,s3c.jiy,x) for 3 < j < 6 . 



(9.3.2) 

(9.3.3) 
(9.3.4) 

(9.3.5) 
(9.3.6) 

(9.3.7) 



(Of course, Xt.s3c.i = 1 is trivially symmetric under x <-!■ y.) The XT.s3c.j for j = 7, 8, 9 are the roots of the 
cubic equation 

- {x + l){x + y + l)e^ + xy{x^ + 2x + y + 2)^ - y^x^ = . (9.3.8) 
The At,s3cj for 10 < j < 15 are roots of the sixth degree equation 



(9.3.9) 
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where 



/eiT = ~{xy + 5 + 5x + 2x^ + Ay + y^) 



(9.3.10) 



/62T = 2a;2y^ + 2x^y + 2 + Ty^x + Ibyx + y^x + lOya;^ + 12a;^ + 10a; + x'^ 

+6x^ + V + + 5y (9.3.11) 

/ggr =3 -~x{2y^x^ + y^x^ + 20?;x^ + 9y + 15j/^a; + lOx^y^ + 4?;^a; + ITyx 

+a;''y + Ix^y + 3 + 5a; + Sx^ + + + 5x^ + lly^ + 6y^) (9.3.12) 

/64T = yx^(5y + 8a; + 2 + 14ya; + 9a;^ + y^x^ + 2y^a; + 4a;^y^ + 4y^ + y^ + a;'' 

+5a;^ + lOya;^ + Oy^a; + 3x^y) (9.3.13) 

/ggj. = -a;^y^(3 + 4y + 2y'^ + y^x + yx^ + + 3x + 3ya;) (9.3.14) 

/66T = x^y^ ■ (9.3.15) 

The next term is 

At,s3c,i6 = xy . (9.3.16) 
Finally, the At.sScj" for 17 < j < 20 are roots of the quartic equation 

+ hire + /42Tf + /43T^ + /44T = (9.3.17) 

where 

/41T = -(3 + yx + 4a; + 2x^ + 3y + a;^ + y^) (9.3.18) 

/42T = a;(2y^x + 4a; + 3 + 4y^ + 3ya;^ + 4yx + 5y + a;^ + 3a;^ + a;^y^ + x^y + y^) (9.3.19) 

/43y = -yx^{2x + 1 + 2ya; + 2ya;2 + 2y + y^ + x'^y^ + x^ + Sa;^) (9.3.20) 

/44T = y^x^ ■ (9.3.21) 
The algebraic equations that yield the twenty terms can thus be summarized as 

{3(1), 2(2), 1(3), 1(4), 1(6)} for sq,3 x L,, FBCy, {T)PBC, (9.3.22) 

by which we mean that these equations consist of three linear, two quadratic, one cubic, one quartic, and 
one sixth-degree equation. This is compared with the results for the Potts model partition functions for the 
corresponding strips of the square lattice with Ly — 1 and 2 in Table ^. 
It is convenient to extract a common factor from the coefficients: 

CT,G,, = . (9.3.23) 

a; — 1 

Of course, although the individual terms contributing to the Tuttc polynomial arc thus rational functions of 
X rather than polynomials in a;, the full Tutte polynomial is a polynomial in both x and y. We have 

CT,Gj = cz,G,j V J . (9.3.24) 
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Table 4: Some properties of Z(G,q,v) for strips of the square lattice of width Ly and arbitrary length, with periodic 
longitudinal boundary conditions. 



Z or P 


Ly 


BCy 




Nx 


eqs. 


rcf. 


ZiG,q,v) 


1 


FBCy 


{T)PBC^ 


2 


{2(1)} 


- 




2 


FBCy 


{T)PBC^ 


6 


{2(1),2(2)} 


[14] 




3 


FBCy 


{T)PBC, 


20 


{3(1),2(2),1(3),1(4),1(6)} 


here 




2 


PBCy 


iT)PBC^ 


6 


{2(1),2(2)} 


here 




3 


PBCy 


PBC^ 


20 


{4(1),3(2),2(3),1(4)} 


here 




3 


PBCy 


TPBCa, 


12 


{3(1),1(2),1(3),1(4)} 


here 


P{G,q) 


1 


FBCy 


{T)PBC^ 


2 


{2(1)} 






2 


FBCy 


{T)PBC^ 


4 


{4(1)} 


[25,24] 




3 


FBCy 


{T)PBC^ 


10 


{5(1),1(2),1(3)} 


[46,47,50] 




4 


FBCy 


{T)PBC^ 


26 


{4(1),1(2),2(3),1(4),2(5)} 


[54] 




3 


PBCy 


PBC, 


8 


{8(1)} 


[51] 




3 


PBCy 


TPBC^ 


5 


{5(1)} 


[51] 




4 


PBCy 


PBC^ 


33 


{9(1),6(2),4(3)} 


[56] 




4 


PBCy 


TPBC^ 


22 


{7(1),3(2),3(3)} 


[56] 



Our calculations are in accord with the following generalization for the family of strips denoted by 

G = sq, Ly X L,,FBCy, {T)PBC,: 

Xt.,g,i = 1 (9.3.25) 

Xz,G,2 - X (9.3.26) 

with coefficients for the cyclic strips 

CT,G.,..,1 =CZ,G.,..,1 ^c^^"^ (9.3.27) 

CT,G.,.. ,2 = cz,G.,.. ,2 - c^^"-!) . (9.3.28) 
For the Mobius strips, the coefficients are determined from the general formulas given in [p5| . 

9.4 Ly = 2 Strips of the Square Lattice with Torus and Klein Bottle Boundary 
Conditions 

We calculate the Tutte polynomials 

6 

T{sq,2x m,PBCy,PBC^,x,y)=Y.''T..-^t,j{XT,s2t.j{x,y)r (9-4.1) 

and 

6 

T{sq,2x m,PBCy,TPBC^,x,y) =Y,CT..s2k,j{\T..s2t^j{x,y)r (9-4.2) 

where 

XT,s2t.i - 1 (9.4.3) 

XT,,2t,2 = X (9.4.4) 
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A 



T,s2t,(3,4) 



2 + 2y + y2 ± 



where 



i?T234 = (x + 2 + 2y + 2/2)2 _ 4^^2 



A 



1 



T,s2t,(5,6) 



where 



Note that 



-fi'T256 = (y^ + 2/ + 1 + 3xy + X + a;2)(2/2 _|_ y _|_ i _ -|_ a; -|_ ■) _ 



AT,s2t J (a;, y) = AT,s2t j (y, a;) for j = 5, 6 . 
The corresponding reduced coefficients (defined as in eq. ( ^.3.23| )) are given by ( 3.3.f(]| )-( 3.3.f6| ). 



(9.4.5) 

(9.4.6) 
(9.4.7) 

(9.4.8) 
(9.4.9) 



9.5 Ly = 3 Strips of the Square Lattice with Torus and Klein Bottle Boundary 
Conditions 



We calculate the Tutte polynomials 
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T{sq, 3 X TO, 

PBCy,PBC^,x,y)=Y, CT.s3t,j{>^T,s3t,j{x, y)y 



and 



T{sq, 3 X m, PBCy,TPBC^,x, y) - ^ CT,s3kj (At ,s3k,j ('^j 
where (ordering the Ar.sstj's by decreasing degree of the associated coefficients cr^sstj] 



AT,s3tj = At,s3j for 1 < j <6 



A 



f 



Z,s3t,{7,8) 



and 



x + y + 4:± V(a^ry + 4j2^^^ 

AT,s3t,9 — y 

AT,s3t,io = xy 



AT,s3tj = AT,s3,i-4 for j = If, f2, f3 



The Xr.sstj for f 4 < j < f 7 are roots of the quartic 



where 



^ + 541tC + .942T^ + 543T^ + 544T = 

541T = -(7 + yx + y^ + x'^ + 5x + 7y + 3y'^) 
542T = y{8x + 6?; + 6 + 4t/2 + + y^ + y^x + iyx^ + ^x^ + 7yx + by'^x + y'^x^) 
g43j. = -y'^x{ix + a;2 + iyx + i + &y + by"^ + y^ + 2y^x + y^x^) 

g44T = y^x^ ■ 



(9.5.f) 
(9.5.2) 

(9.5.3) 

(9.5.4) 

(9.5.5) 
(9.5.6) 

(9.5.7) 

(9.5.8) 

(9.5.9) 
(9.5.f0) 
(9.5.ff) 
(9.5.f2) 
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Finally, the Ar.sst.j for j = 18, 19, 20 are roots of the cubic equation 



where 



331T 



C + ffSirC + + 533T — 



4 + 5(x + y) + 3(a;^ + y^) + 2a;?; + + 



ff32T = xy 



1 + 2(a; + y) + 3(x^ + xy{x + y) + y^) + Axy + (x'^ + y'^) + x^y^ 



533T 



4 4 

-y X 



Among the Ar.sst.j's that depend on x and y, ten are symmetric functions of these variables: 
>'T,s3tj{x,y) ^ XT,s3tj{y,x) for j = 3,4,5,6, 7,8, 10,18, 19,20 . 



(9.5.13) 

(9.5.14) 

(9.5.15) 
(9.5.16) 

(9.5.17) 



(and Xx.sst,! = 1 is trivially symmetric under x <-> y.) The terms that enter in eq. (9.5.2) for the strip with 
Klein bottle boundary conditions are 

AT,s3fcj = AT,s3t,j for j = 1,2 (9.5.18) 

XT,s3k.j = AT,s3tj+4 for 3 < j < 5 (9.5.19) 

AT,.3fej = AT,.3t,j+8 for 6 < i < 12 . (9.5.20) 

The eight terms in the Tutte polynomial for the torus strip \T,s3t.j with 3 < j < 6 and 10 < j < 13 do not 
occur in the Tutte polynomial for the Klein bottle strip. 

For both of the types of strips = sg, 3 x m, PBCy, {T)PBCx, the corresponding coefficients satisfy 

ct,Gs,] = cz.c.j ^J- (9.5.21) 

9.6 Special Values of Tutte Polynomials for Ly = 3 Strips of the Square Lattice 
9.6.1 General Relations 

For a given graph G ~ {V,E), at certain special values of the arguments x and y, the Tutte polynomial 
r(G, X, y) yields quantities of basic graph-theoretic interest. We recall some definitions: a spanning subgraph 
was defined at the beginning of the paper; a tree is a connected graph with no cycles; a forest is a graph 
containing one or more trees; and a spanning tree is a spanning subgraph that is a tree. We recall that the 
graphs G that we consider are connected. Then the number of spanning trees of G, Nst{G), is 



7VsT(G) = r(G,l,l) , 
the number of spanning forests of G, Nsf{G), is 

7V5F(G)=r(G,2,l) , 
the number of connected spanning subgraphs of G, NcssciG), is 

NcssG{G)=TiG,l,2) , 
and the number of spanning subgraphs of G, Nssg{G), is 

iVssG(G) =T(G,2,2) . 



(9.6.1) 



(9.6.2) 



(9.6.3) 



(9.6.4) 
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An elementary theorem (e.g., [|T6|) is that 

Nssg{G) = 2^(^) . (9.6.5) 

Since T{G,Ly x m,x,y) grows exponentially as m ^ oo for lattice strips G of the type considered here, 
for the values {x,y) — (1, 1), (2,1), (1,2), and (2,2), one defines the corresponding constants 

Zset{{G})= lim n{G)-'^ In NsetiG) , set = ST, SF, CSSG, SSG (9.6.6) 

n(G)^oo 

where, as above, the symbol {G} denotes the limit of the graph family G as n(G) oo. These constants 
are determined completely by the term Xz,G,j that is dominant in the PM phase. Since this is the same 
for a given G with a specified transverse boundary conditions, independent of the longitudinal boundary 
conditions, the resultant value of Zset is the same for a strip of a given lattice type A, width Ly, and transverse 
boundary condition, independent of the longitudinal boundary condition {FBC'x, PBCx,TPBCx) @. 

9.6.2 Free Transverse Boundary Conditions 

For the Ly = 3 cyclic and Mobius strips of the square lattice we thus have 



zsT,3xoo.FBCy = ^ ln[AT,.3c,i7,.=i,a=i] = 0.84307429... (9.6.7) 



As a special case of the result noted above, this is the same for the open Ly — 3 strip of the square lattice. 
A general upper bound on the number of spanning trees for a graph G is 



1 / 2\E\ ^ ""^ 



NsTiG)<-[^^] . (9.6.8) 
n \n — 1 



For n oo, in terms of the effective degree A^fj defined in ( 1.29 ), this yields 

2ST,{G} < In A,/^ . (9.6.9) 
That is, for the cyclic, Mobius, or free strip of the square lattice with width Ly, 



ZST,sq,LyXoo,FBCy, < In 



4fl-— 1 

V 2LyJ 



(9.6.10) 



For Ly = 3, this evaluates to 



zsT,sq,3x^,FBCy < ln(10/3) = 1.2039728... (9.6.11) 
The ratio of the actual zst to this upper bound is thus 

= 0.7002436... (9.6.12) 



ZST,sq,3xcx.FBCy 



ln(10/3) 



We may also compare the value zst in ( 3.6.7 ) with the value for the Ly — 2 cyclic, Mobius, or free strips of 
the square lattice, namely (e.g., p^) 

ZST,sq,2xoc,FBCy = - ln(2 + \/3) 0.6584789... (9.6.13) 
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Thus, as one increases the width of the strip with FBCy from Ly = 2 to Ly = 3, the value of zst increases 
about 30 %. Another comparison of interest is the ratio of zst for these finite- width strips with zst for the 
fuU 2D square lattice, which has the value Q 



ZST.sg = -Gcat. = 1.1662436. 



where Gcat. is the Catalan constant, 



Gcat. = ^ 



{-ly 



(2j + 1)^ 



We have 



and 



ZST,sq.2y.oa,FBCy 
ZST,sq 



Zst, sn.Sxoo, FBCy 



= 0.91596559. 



0.564615... 



0.722897. 



ZST,sq 



(9.6.14) 



(9.6.15) 



(9.6.16) 



(9.6.17) 



Thus, for Ly ~ 3, the value of zsT,sq.LyXoo,FBCy is within about 30 % of the value (9.6.14) for the infinite 
square lattice. 

Our calculations also yield, for spanning forests. 



ZSF,3xoc,FBCy = ^ ln[AT,s3c,17,a;=2,y=l] = 1.106676. 



and, for the number of connected spanning subgraphs, 

1 



ZCSSG,3xcc,FBCy = -^M^T,s3c,17,x=l,y=2\ = 0.966005. 



From the general theorem (9.6.5), together with the definition (1.28), we have 



ZSSG 



ln2 . 



(9.6.18) 



(9.6.19) 



(9.6.20) 



(9.6.21) 



In the present case, with Ag// — 10/3 for the Ly = 3 cyclic or Mobius strip graph, we have 

^ssG = ^ln2 = 1.155245... 
9.6.3 Periodic Transverse Boundary Conditions 

For the Ly — 2 and Ly = 3 strips of the square lattice with periodic transverse boundary conditions (torus, 
Klein bottle and cylindrical b.c.) we have 



zsT,2xoo,PBCy = ^ ln[AT,s2t,5,2:=i,y=i] = ^ ln(3 + 2\/2 ) = 0.8813735. 



(9.6.22) 



ZST,3xoo,PBCy — 



— ln[AT,s3i,18,x=l,j/=l] ~ 2 ^-"^ 



23 + 5V21 



In 



= 1.0445328. 



(9.6.23) 
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These results agree with the eq. (6.6.1) in |Q, where a general discussion was given for the number of 
spanning trees on various lattices and lattice strips. The torus is a A-regular graph and, even for the 
cylindrical strip, in the ^ oo limit, the effect of the end vertices goes to zero, so that in all cases, the 
coordination number is 4. Hence, the upper bound (9.6.8) reads 

ZST.sq,LyXoo,PBCy < In 4 . (9.6.24) 

The ratios of the actual values of zst for Ly = 2 and = 3 to this upper bound are thus 

zsT,,q.2xoo,PBCy ^ 0.635777... (9.6.25) 
ln4 

zsT,sq,3xoo.PBCy ^ 75347^^^^ g 26) 

In 4 

Although the Ly = 3 strips with free transverse boundary conditions are not A-regular graphs, the cor- 
responding strips with periodic transverse boundary conditions (and, say, periodic longitudinal boundary 
conditions) are A-regular, so one may also compare the exact results with a sharper upper bound for a 
A-regular graph G with n vertices [|o[ |l| , 

^STiG) < ^-^yC.r (9.6.27) 

where 

(A - 1)^-1 
- [A(A-2W-i 

which yields 

2{G} < ^{G},A =lnCA (9.6.29) 
where we label the upper bound for the limit of A-regular graphs as z^Qy We have 

ZST,2xoo,PBCy 



31n(3/2) 

ZST,3xoo,PBCy 

31n(3/2) 



0.724578... (9.6.30) 
.8587116... (9.6.31) 



Comparing with the value of zst for the 2D square lattice (cf. eq. ( 9.6.14 ) above), we have 

ZST,sq.2xoo,PBC. 



ZST,sq 
ZST,sq,3xoo,PBCy 



= 0.755737... (9.6.32) 
= 0.8956386... (9.6.33) 

ZST,sq 

It is impressive that the strips with rather modest Ly = 3 width and periodic transverse boundary conditions 
(and any longitudinal boundary conditions) yield a value of zsT,sq,LyXoo,PBCy that is only about 10 % below 
the value for the infinite 2D square lattice. 
We also calculate 



ZSF,2xoo,PBCy — 2 l^['^T,s2t,5,a;=2,y=l] — 2 ^^^^ 



i(ii + VT05) 



= 1.181533... (9.6.34) 



ZSF,3xoo,PBCy = X ln[AT,s3t,18,a:=2,2/=l] = 1.263294... (9.6.35) 
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ZCSSG,2xoo,PBCy — 2 l^['^T,s2t,5,2; = l,y=2] = '^^l-'^ 



-(11 + a/105) 



= 1.181533... 



and 



ZCSSG,3xoa,PBCy — g ln[AT,s3t,18,2; = l,a=2] 



1.263294.. 



(9.6.36) 



(9.6.37) 



The equality of ( |j.6.36 ) with (9.6.34) and the equaUty of (9.6.37) with (9.6.35) follow from the x y 
symmetry of XT.s2t,b a-nd Xt.sM,!?,, as indicated above in eqs. (|X|) and ( 19.5.171 ). Finally, from eq. ( ^X20| ) 
with A = 4, we have, for both the torus or Klein bottle strip with arbitrary Ly, 



zssG = 21n2 = 1.386294... 



(9.6.38) 
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